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Abstract We review the properties and nature of luminous high-redshift ra- 
i-C ■ dio galaxies (HzRGs, z > 2) and the environments in which they are located. 

HzRGs have several distinct constituents which interact with each other - 
relativistic plasma, gas in various forms, dust, stars and an active galactic 
nucleus ( AGN) . These building blocks provide unique diagnostics about con- 
ditions in the early Universe. We discuss the properties of each constituent. 
Evidence is presented that HzRGs are massive forming galaxies and the pro- 
genitors of brightest cluster galaxies in the local Universe. HzRGs are located 
in overdense regions in the early Universe and are frequently surrounded by 
protoclusters. We review the properties and nature of these radio-selected 
protoclusters. Finally we consider the potential for future progress in the 
field during the next few decades. A compendium of known HzRGs is given 
| in an appendix. 

Keywords Radio galaxies • High-redshift • Massive galaxies ■ Clusters 

o 

OO 

O 



G. K. Miley 

Sterrewacht, Leiden University, Postbus 9513, 2300RA Leiden, The Netherlands 

Phone: +31 71 5275849 

E-mail: miley@strw.leidenuniv.nl 

C. De Breuck 

European Southern Observatory, Karl-Schwarzschild-Strasse 2 D-85748, Garching, 
Germany 

Phone: +49 89 32006613 
E-mail: cdebreuc@eso.org 



2 



George Miley, Carlos De Breuck 



Contents 



1 INTRODUCTION 

1.1 Scope of this article - HzRGs 

1.2 History 

1.3 Hunting for HzRGs 

1.4 Redshift distribution of HzRGs 

1.5 Constituents of HzRGs 

2 RELATIVISTIC PLASMA 

2.1 HzRGS and low-redshift radio galaxies 

2.2 Radio sizes and morphologies. Size as an evolutionary clock .... 

2.3 Radio size vs redshift correlation 

2.4 Radio spectral index vs redshift correlation 

2.5 Nonthermal X-ray emission 

3 GAS AND DUST 

3.1 Hot ionised gas - Radio depolarization 

3.2 Warm ionized gas 

3.2.1 Ionisation 

3.2.2 Abundances and star formation 

3.2.3 Inner Kinematics: Outflows, jet interactions and superwinds 

3.2.4 Outer kinematics: Infall of the quiescent halos 

3.2.5 Relation to non-radio Lya nebulae 

3.2.6 Nebular continuum 

3.3 Neutral Gas 

3.3.1 HI Absorption 

3.3.2 Lya absorption 

3.4 Molecular Gas 

3.5 Dust 

3.5.1 Millimeter and sub-millimeter emission 

3.5.2 UV continuum polarisation 

4 STARS 

4.1 Old stars: HzRGs as the most massive galaxies 

4.2 Young stars - UV absorption lines 

4.3 The alignment effect - Jet-induced star formation 

5 ACTIVE GALACTIC NUCLEUS (AGN) 

5.1 Hidden or dormant quasar - The Unification Church 

5.2 Supermassive black hole - Powerhouse of the AGN 

5.3 Extinction of the "dinosaurs" 

6 NATURE OF HzRG HOSTS - MASSIVE FORMING GALAXIES . . . 
6.1 The Spiderweb Galaxy - a case study 

7 ENVIRONMENT OF HzRGS - PROTOCLUSTERS 

7.1 Finding protoclusters using HzRGs 

7.2 Properties of radio-selected protoclusters 

7.3 Are radio-selected protoclusters typical? 

7.4 Protocluster evolution 

8 CONCLUSIONS AND FUTURE 

9 APPENDIX: KNOWN HZRGs 



Distant Radio Galaxies and their Environments 



3 



1 INTRODUCTION 

Martin Harwit. " Yes. I think that if there were less requirement for theoret- 
ical justification for building a instrument that is very powerful, you would 
then have a better chance of making a discovery. " 

Bernard Burke. "The record of radio astronomy is clear. Quasars were not 
found by the desire to find black holes or by the need to find long-distance 
cosmological probes. " 

Proceedings of Greenbank Workshop on Serendipitous Discovery in Radio 
Astronomy, 1983 

Distant radio galaxies are among the largest, most luminous, most massive 
and most beautiful objects in the Universe. They are energetic sources of ra- 
diation throughout most of the electromagnetic spectrum. The radio sources 
are believed to be powered by accretion of matter onto supermassivc black 
holes in the nuclei of their host galaxies. Not only are distant radio galaxies 
fascinating objects in their own right, but they also have several properties 
that make them unique probes of the early Universe. 



1.1 Scope of this article - HzRGs 

Our review will be restricted to galaxies that have redshifts, z > 2 and 
radio luminosities at 500 MHz (rest frame) L 50 o( r est) > 10 27 5 W Hz -1 . Radio 
emission from such objects has a steep nonthermal spectrum, is collimated 
and is usually extended by tens of kiloparsec. There are very few known z > 
2 radio galaxies associated with compact flat-spectrum radio sources. 

This definition is somewhat arbitrary. Luminous distant radio galaxies 
have properties that are different from less powerful radio galaxies at low red- 
shifts and the properties of radio galaxies change gradually with luminosity 
and redshift. Likewise, z > 2 quasars that are associated with steep-spectrum 
extended radio emission have many similar properties to our distant radio 
galaxies (see Section [5. II) and are located in similar environments. However, 
without some restriction, our task would not have been tractable. 

For concisenes s, we shall f r equen tly refer to such high-redshift radio galax- 
ies as "HzRGs". iMcCarthyl (|l993l ) gave an extensive review of this topic 
in 1993. Since then there have been several workshops on high-redshift ra- 
dio galaxies, whose pro ceedings have been published, includi ng Amsterdam 
( Rotteering et all 1999h. L eiden ( Jarvis and Rottgeringl [20031 ) and Granada 



(|Villar-Martm et all 120061 ). Thousands of refereed papers dealing with dis- 
tant radio galaxies and related topics have been published during the last 
fifteen years. Although our goal is to be comprehensive, some personal bias 
will inevitably have entered in selecting the topics and literature. Our apolo- 
gies for this. 

The structure of this review is as follows. We shall first set the scene by 
reviewing the history of the field. Then we shall describe techniques used for 
finding HzRGs and discuss how HzRGs are distributed in redshift. Following 
an overview of the diverse constituents of HzRGs, we shall discuss each HzRG 
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component in detail. The nature of HzRGs and their role in the general evo- 
lution of galaxies will then be covered. We shall present evidence that HzRGs 
are the massive progenitors of dominant cluster galaxies in the local Universe 
and that they are located in forming galaxy clusters. The properties of these 
radio-selected protoclusters are then reviewed. Finally, we shall discuss some 
of the exciting prospects for future HzRG research. An appendix includes a 
list of known HzRGs at the time of writing (October 2007) . 



1.2 History 

The study of distant radio galaxies has progressed in several phases: 
(i) Infancy. Mid-1940s to mid-1960s: After the discovery that Cyg nus A 



(ij Infancy. Mid-iy^Us to mid-Wbus: Atter the discovery that (Jyg nus A 
is associated with a faint distant galaxy ( Baade and Minkowskill 19541 ), radio 



astronomy became one of the most important tools of observational cos- 
mology. Redshifts of up to z = 0.45 were measured for galaxies associated 
with radio sources ( Minkowskilll960f ) . Two fundamental breakthroughs were 



made towards the end of this period. First, quasi-stellar radio sources or 
" quasars" were discovered, with much larger redshifts even than radio galax- 
ies. Secondly, it was demonstrated that the space density of radio sources 
varies with cosmic epoch, an observation that sounded the death knell for 
the Steady State cosmology 

(ii) Childhood. Mid-1960s to mid-1980s: This can be termed the "Spinrad 
era" . Hy Spinrad devoted a considerable amount of his time, and that of the 
120-inch telescope at Lick, to measuring the redshifts of faint radio galaxies 
by means of long photographic exposures. With c onsiderable effort, he pushed 
the highest radio galaxy redsh ift out to z ~ 1 ( Spinradlll97rj : ISpinrad et all 



119771 : ISmith and Spinradl fl980). At that time, it was not realised that radio 
sources interact with their optical hosts and the different wavelength regimes 
were usually studied in isolation from one another. Radio sources were re- 
garded as interesting high-energy exotica that were useful for pinpointing 
distant elliptical galaxies. They were not believed to play an important role 
in the general scheme of galaxy evolution. During this period the paradigm 
that radio sources are powered by accretion of matter onto rotating black 
holes was developed. 

(Hi) Teens. Mid-1980s to mid-1990s: The replacement of photographic tech- 
niques by CCDs revolutionised optical astronomy. Large numbers of very 
distant radio galaxies were discovered with redshifts out to z ~ 5 (HzRGs). 
The surprising discovery that the radio sources and the optical host galaxies 
are aligned showed that there must be considerable interaction between the 
radio sources and their host galaxies. The role that orientation could have in 
determining observed properties was realised." Orientation unification" that 
explained differences in the observed types of active galaxies as merely due to 
the viewing angle of the observer became a popular interpretative " do-it-all" . 
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(iv) Maturity. Mid-1990s to the present: New efficient techniques for find- 
ing distant galaxies were developed. These included photometric searches for 
"dropout" objects with redshifted Lyman break features and narrow-band 
searches for objects with excess redshifted Lyman a fluxes. The advent of 
these new techniques placed the cosmological use of radio galaxies in a new 
perspective. Although radio astronomy lost its prime position as a technique 
for finding the most distant galaxies, it became clear that luminous radio 
galaxies play a highly important role in the evolution of galaxies and the 
emergence of large scale structure. The discovery of a relation between the 
masses of elliptical galaxies and the inferred masses of nuclear black holes 
led to the conclusion that all galaxies may have undergone nuclear activity 
at some time in their histories. This brought nuclear activity into the main- 
stream of galaxy evolution studies. 



1.3 Hunting for HzRGs 

Appendix A contains a compendium of 178 radio galaxies with z > 2 known 
to the authors. This list is the result of more than two decades of hard work. 
The relatively small number of known HzRGs illustrates both their rarity 
and the difficulty of detecting them. 

Finding distant radio galaxies involves a multistage process. 

— The first st ep (e.g. lRottgering et allll994l ; lBlundell et aJl998HDe Breuck eta] 
200(J l2004J l is to filter out probable HzRGs from the huge number of ra- 
dio s ources contained in low-frequency radio surveys (e.g. ICondon et all 
1998). The usual criteria used for filtering likely candidates are (i) radio 
colours - extremely steep spectra (Section l2.4j) and (ii) small angular sizes 
(Section CO]). 

— In the second step, relatively bright nearby objects are discarded, by 
comparing the radio positions with existing wide-field shallow surveys at 
optical and infrared wavelengths. 

— The third step is to refine the positions of the remaining radio sources 
to arc-second accuracy, thereby facilitating their identification with faint 
galaxies. Until t he 1990s, radio identifications were made u sing optical 
CCD data (e.g. iRottgering et all Il995k iMcCarthv et~aJll996l ). However, 
after sensitive near-IR detectors became available, the near-IR K— band 



was found to be more efficient for identifyinj 


? distant radio sources (e.g. 


Jarvis et alll200lHDe Breuck et all2002ll2004; 


Jarvis et alll2004l:ICruz et all 


2006). More than 94% of bright radio sources are identified down to 



if =22. Furthermore, the if— band magnitude of an HzRG provides a 
first indication of its redshift, by means of the K — z relationship (see 
Section l4~Tj) . 

— The fourth step is to carry out spectroscopic observations of the HzRG 
candidates with large optical and infrared telescopes. These are needed 
to determine redshifts from the wavelengths of HzRG emission lines (e.g. 
Rottgering et a]|l997tlpe Breuck et all200ll;|Jarvis et afeOQlHDe Breuck et all 
2006t iBornancini et a]|l2007[ ) (see Section l3~2|) . 
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The main limitation on finding HzRGs has been the scarcity of available 
observing time on large optical/IR telescopes for the last stage in the process. 
Because candidates are located too far apart in the sky to permit multi- 
object spectroscopy, tedious long-slit spectroscopy on many individual fields 
is required to determine the redshifts. 

Although radio selection ensures that there is no a priori selection against 
dust properties (e.g. extinction), there are other observational selection ef- 
fects that introduce bias into the redshift distributions of HzRGs and the list 
given in the appendix. The first obvious source of bias is that the determina- 
tion of a spectroscopic redshift is dependent on being able to observe bright 
emission lines with gro und-based telesco pes. The primary line for such red- 
shift measurement (e.g. lMcCarthvlll993| ) is Lya A1216A (typical equivalent 
widths of several hundred A). Useful additional lines are CIV A1549A, Hell 
A1640A and CIII] A1909A (equivalent widths in excess of ~ 60 A). There is a 
so-called "redshift desert" 1.2 < z < 1.8, fe.g. lCruz et"alll2006l )h where [Oil] 
A3727A is too red and Lya A1216A is too blue to be easily observed from 
ground-based optical telescopes. Objects located within the redshift desert 
are thus under-represented in samples of HzRGs. 

A second source of bias is redshift incompleteness. A small fraction of 
radio sources, (~ 4%) are not identified to K ~ 22 and about a third of 
those with K-band identifications do not show any emission or absorption 
lines, even after long exposures on 8 - 10m class optical/IR telescopes. These 
objects either (i) have emission lines that are hidden by substantial dust 
obscuration (e.g. |Pe Breuck et a i l200ll:lReuland et all [2003). (ii) are located 
at such large redshifts that Lya and other bright emission lines fall outside the 
easily observable spectral windows, or (iii) are peculiar in that they radiate 
no strong emission lines. 



1.4 Redshift distribution of HzRGs 

HzRGs are extremely rare. Luminous steep-spectrum radio sources associ- 
ated with HzRGs have typical luminosities of \ji.ighz > 10 33 erg s^ 1 Hz -1 
ster -1 . The number density of radio sources with this luminosity in the red- 
shift range 2 < z < 5 is a few times 10~ 8 Mpc~ 3 , with large uncer tainties 
(|Dunlop and Peacock|[T990t IWillott et alll200lt IVenemans etall2007| ). 

Although HzRGs are sparsely distributed in the early Universe, such ob- 
jects are almost nonexistent at low redshifts. The co-moving space density of 
luminous steep-spectrum radio sources increases dramatically by a factor of 
100 - 1000 between < z < 2.5 and then appears to flatten out( Willott et all 



l200lfc IJarvis et alll200lh . Although there are huge uncertainties in the evolu- 
tion of the luminosity function at higher redshifts, no significant cut-off in 
space density has yet been observed. 

The redshifts at which the space density of radio sources is maximum 
correspond to a crucial era in the evolution of the Universe. It coincides with 
the epoch when (i ) luminous quasars also ap pear to have had their maximum 
space density (e.g |Peilll995t iFan et alll200lf ) (ii) star formation was rampant 
and more than an order of magnitude larger than the present (jMadau et all 
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1998: iLillv et alllT996t ISchiminovich et afeOOSh and (iii) galaxy clusters were 
forming, but were not yet gravitationally bound structures. 

The question of why the most luminous radio sources (and the most 
luminous quasars) became extinct in the local Universe is an intriguing one 
that is still not understood (Section 15. 3j) . 

1.5 Constituents of HzRGs 

Radio galaxies have several distinct emitting components which provide di- 
agnostics about various physical constituents of the early Universe. A list of 
known HzRG building blocks is given in Table 1 , together with a summary 
of techniques used to study them. Also included are a list of the resultant 
diagnostics, some useful relevant references and our best estimate for the 
typical mass of the component in HzRGs. We caution the reader to thor- 
oughly examine the assumptions inherent in deriving each of the diagnostics. 
Although, these assumptions are needed in order to reach some conclusions 
about the nature of HzRGs, it is important never to forget that they exist 
and to be sceptical about the "houses of cards" that are frequently involved. 

Because they are highly luminous and (unlike quasars) spatially resolvable 
from the ground, most components of HzRGs provide important diagnostic 
information about the spatial distributions of processes within HzRGs and 
their environment. The fact that the different constituents are present in the 
same objects and that the interrelationships and interactions between them 
can be studied make distant radio galaxies unique laboratories for probing 
the early Universe. 

As can be seen in Table 1, several constituents of HzRGs are inferred to 
be extremely massive, including old stars (up to ~ 10 12 M Q ), hot gas (up to 
~ 10 12 M ) and molecular gas (up to ~ 10 11 M Q ). 

Figure [1] shows the spectral energy distribution (SED) of a typical HzRG 
from radio to X-ray wavelengths, together with a decomposition into various 
observable HzRG constituents - relativistic plasma, gas and dust, stars and 
the active galactic nuclei (AGN). We shall discuss each of these building 
blocks individually in Sections to Note from Figure [1] that disentangling 
the various components in the optical and the infrared is difficult and the 
results are often extremely model dependent and several alternative solutions 
may be equally consistent with the available multi-wavelength data. 

2 RELATIVISTIC PLASMA 

"All things are double, one against another" Ecclesiasticus (xlvii. 24) 

"The problem is how much to stress the sameness of twins and how much 

to emphasise their differences.", Baby and Child Care, Benjamin Spock, 1946 

Shortly after the discovery of extragalactic radio sources, their nonther- 
mal spectra and polarisations led to the conclusion that their emission is pro- 
duced by synchrotron radiation from a relativistic plasma. The large radio 
luminosities together with the lifetimes of the synchrotron-emitting electrons 
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Table 1 

Constituents of distant radio galaxies 



PrvTl ci" it'll f^iri" 


Ohtiprva r i 1 o 


' 1 Vm f*n 1 Tl 1 n crn risst 1 i^k; T-?P»fc; 

1 V J 1 1 k 1 1 1 UldgllUSLll^O 1 V < 1 — . 


Mass 

(Ms) 


Kelativistic plasma 


Radio continuum 


Magnetic tield, age, en- 1,2 
ergetics, pressure, par- 
ticle acceleration. Jet 
collimation and propa- 
gation 




X-ray continuum 


Magnetic tield, equipar- 3,4,1 
tition, pressures 


Hot ionized gas 
T e ~ 10 7 -10 8 K 
n e ~ 10 _1 ' 5 cm~ 3 


Radio 

(de) polarisation 


Density, magnetic lield, 1 


10 u-i2 


X-rays 


Temperature, density 
mass 


Warm ionized gas 
T e ~ 10 4 -10 5 K 
n e ~ 10°- 5_1 - 5 cm~ 3 


UV-optical 
emission lines 


Temperature, density, 5,6,7,8 
kinematics, mass, ion- 
isation, metallicity, 
filling factor 


^q9-10.5 


Nebular 
continuum 


SED contamination 9,10 


Cool atomic gas 
T s ~10 3 K 
n(HI) ~ lO^m" 3 


HI absorption 


Kinematics, column 11,8 
densities, spin temper- 
ature, sizes, mass 


10 7-8 


U V-optical 
absorption lines 


Kinematics, mass, col- 8,12 
umn densities, metallic- 13,14 
ity 


Molecular gas 
T ~ 50 - 500K 
n(H 2 ) > 10 2 cm" 3 


(Sub)millimeter 
lines 


Temperature, density, 15 
mass 


10 io-n 


Dust 

T - 50 - 500K 


UV-optical 
polarisation 


Dust composition, scat- 16 
tering, mass, hidden 17 
quasar 


10 8-9 


(Sub)millimeter 
continuum 


Temperature, mass, 18 
heating source 


Old stars 
t > 1 Gyr 


Optical to near 
IR 

continuum 


Age, mass, formation 19 
epoch 


10 U-12 


Young stars 
t <0.5Gyr 


UV-optical 


Star formation rates, 20,8 
ages, history 


10 9-10 


Lya 


Star formation rate 20 


Quasar (hidden or 
dormant) 


UV-optical 
polarisation 
broad lines 


Luminosity 21,22 




Supermassive black 
hole(SMBH) 


Extended radio, 
Quasar 


Formation, evolution 23,24 


~10 9 



□Refer ences : 1 = IMileyl lfl980T) . 2 =IKlamer et all (I2006T). 3 =IFelten and Morrison 
11966T) . 4 = ISchwartzl (120021) . 5 =lQsterbrock and Ferlandl (12009 ) 6 = IGroves et a 
■■>"»*, 7 = IGroves et all <I2004I). 8 = IDopita and Sutherland! 1(2003). 9= lAller 
, 10 = IDickson et all (119951 '), 11 = iMorgantil (I200a)7 l2 = Ivan Oiik et a 
, 13 = Bi nette et all (120001) . 14 = IBinette et all (120061. 15 = IDownes eta 
, 16 = Cimatti et aj (THm) . 17 = IVernet et d (l2fJ0lT) 18 = IReuland eta 
, 19 = Sey mour o) al (120071) . 20 = IMadau et all TSU ), 21 =TCimatti et a 
, 22 = IVernet et all l|2001l ). 23 = IBlandfordl l|2001l ). 24 = IPeacockl l| 19991 ) 
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Restframe Wavelength (/jm) 
10- 4 10- 3 0.0 1 0.1 1 10 1 00 1 00 1 4 1 5 1 6 




Observed Wavelength (/i.m) 

Fig. 1 Spectral energy distribution (SED) of the continuum emission from the 
HzRG 4C23.56 at z = 2.5, illustrating the contributions from the various con- 
stituents. [From De Breuck et al. in preparation]. Coloured lines show the decom- 
position of the SED into individual components, under many assumptions. Cyan = 
rad io synchrotron (Section[2]). Black = Abso rbed n onthermal X-ray AGN (Section 
12. 5p . Yellow = nebula r con tinuum (Section 13.2.6)1 . Blue = AGN-heated th erma l 
dust emission (Section 13.5) 1 . Red = Starburst- heated dust emission (Section 13,51 1, 
Green = Stars (Section |4|. Magenta = scattered quasar (Section [5]). The addition 
of the overlapping modeled components fits the SED well. 

imply total energies of > 10 60 erg (e.g. lMilevlflQSOh . These huge energies, the 
collimation on sub-parsec scales and the similarity of orientation between 
the compact and overall radio structures led to the conclusion that the colli- 
mated relativistic beams are produced by rotating supermassive black holes 
(SMBHs) located at the centres of the host galaxies. The SMBHs are pos- 
tulated to have rotation axes aligned with the radio source axes. The radio 
sources are powered through gravitational energy from material accreting 
onto the SMBHs, that is converted into kinetic energy of the collimated rel- 
ativistic jets (See Section [5T2]) . 

2.1 HzRGS and low-redshift radio galaxies 

Extragalactic radio sources can be classified according to their sizes. 

— Most radio sources detected in low-frequency radio surveys [y <; 2 GHz) 
are "extended", with projected linear sizes ranging from several tens of 
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kiloparscc to several megaparsecs, i.e. much larger than those of their op- 
tical galaxy hosts. Extended extragalactic radio sources have nonthermal 
spectra, with typical spectral indices at low-redshifts of a ~ -0.7, with 
flux proportional to v a . An example is Cygnus A in Figure [31 

— Most radio sources detected in high-frequency radio surveys are "com- 
pact", with typical sizes of < 1 pc and relatively flat spectral indices (a 
> -0.4). 

— A small fraction of sources have linear sizes between ~ lOOpc to a few 
kpc. These sources have peaked spectra that are self absorbed at low fre- 
quencies (< a few GHz) and are variously known as "peakers" , gigahertz- 
peaked spectrum (GPS) sources or compact steep-spectrum (CSS) sources. 

Galaxies and quasars host extragalactic radio sources of all three classes, 
with extended radio sources predominantly identified with galaxies and com- 
pact sources mainly associated with quasars. 

HzRGs generally host extended extragalactic radio sources and differ in 
several properties from low-redshift radio galaxies. At larger redshifts the 
typical radio luminosities increase, the typical sizes decrease (Section 12.31) 
and the typical radio spectra steepen (Section 12 .4[) . The host galaxies of 
HzRGs differ from the low-redshift radio galaxy hosts in (i) the presence of 
emission-line halos (Section 13. 2\ (ii) increased dumpiness of the continuum 
emission (Section I4.3j) and (iii) alignments of the radio sources with both 
the emission-line halos (Section l3.2.3p and the UV/optical galaxy continuum 
emission (Section I4.3[) . 



2.2 Radio sizes and morphologies. Size as an evolutionary clock 



The radio structure s of H zRGs have been studied by ICarilli et all ( 1997h 



and iPentericci et a 3 (|2000h . Both articles describe observations of 37 ra- 
dio galaxies with z > 2 using the VLA at 4.7 and 8.2 GHz, at resolu- 
tions down to 0.25 arcec. In accordance with their large radio luminosities, 
most HzRGs are " Fanaroff- Riley Class II" radio sources ([Fanaroff and Rileyl 
Il974h . with double structures, edge-brightened mophologi es and one or more 
hot-spots located at the extremit ies of their lobes (|Milevl Il980t ICarilli et all 
11994 Il997t IPentericci et alll2000h . In general HzRGs do not have apprecia- 
ble flat-spectrum core components at their nuclei, but the use of the ultra- 
steep spectrum criterion in searching for HzRGs discriminates against find- 
ing HzRGs w ith flat-spectrum cores. Standard minimum energy assumptions 
(|Milevlll980D give typical minimum pressures in these hotspots of a few x 10 



dyn cm -2 and corresponding magnetic field strengths of a few hundred ^iG. 

There have been a few VLBI observations of fine-sc ale structure in the 
lobes of extended radio sources associated with HzRGs dGurvits et a lll997t 
ICai et all 12001 iPerez-Torres and De Breucld l2005t iPerez-Torres et alll2006£ 
A component of size ~ 65 pc has been detecte d several kpc from the nucleus 
of 4C41.17 at z = 3.8 ( Perez- Torres et a 1 2006). Consideration of the energet- 



ics suggests that the radio component is associated with a gas clump of mass 
Mb > 1.5 x 1O 8 M0. Intriguingly, this is typical for the mas ses of "predis- 
rupted clumps" invoked as the progenitors of globular clusters ([Fall and Reed 
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Fig. 2 Radio structure of 4C41.1 7 at z = 3.8. [From van Breugel and Reuland, 
private communication]. See also (Reuland ct al 2003). Contours obtained at 1.4 
GHz with the VLA are superimposed on a K eck narrow-band image in redshifted 
Lya, showing the warm ionized gas (Section 13.21 The radio angular size of ~ 13 
arcseconds corresponds to a projected linear size of ~ 90 kpc. The radio spectrum 
of this object is shown in Figure[3j;. Figure [5 ] illu strated the radio rotation measure 
of the brightest radio component and FigurefT5]is a higher resolution radio contour 
map of the central region. 



1977). Investigating fine scale radio structure in HzRGs and using the radio 
jet interactions to probe the interstellar medium in the early Universe are 
likely to be an important field of study for the more sensitive long-baseline 
interferometers that are presently under construction, such as e-Merlin and 
the e-EVN. 

There are sev eral correlations betwe en the sizes of radio sources and and 
other properties ( Rottgering et all l2000). First, for 3C sources at lower red- 
shifts, there is a relation with optical morphology. Smaller radio sources con- 
sist of several bright optical knots aligned along th e radio axes, while larger 
sources are less lumpy and aligned ( Best et alll998l) . Secondly, there is a rela- 
tion with emission line properties. Smaller radio sources have generally lower 
ionisatio ns, higher emis sion line fluxes and broader line widths than larger 
sources ([Best et allll999i ) . Thirdly, there is a connection with Lya absorption. 
HzRGs with smaller radio sizes are more likely to h ave Lya absorption than 
larger sources (Section |3.3.2[) ( van Oiik et all fT997[ ). Taken together, these 
correlations are consistent with an evolutionary scenario in which radio size 
can be used as a "clock" that measures the time elapsed since the start of 
the radio activity. 

Note that the large sizes of radio sources (usually several tens of kpc) 
imply that the nucleus of host galaxies has been undergoing activity for 
at least > 10 5 y (light travel time across the source) and up to > 10 8 y 
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(assuming that the jets advance at a few hundred km/s). Hence we know 
that AGN associated with extended radio sources must be long-lived. This is 
not the case for radio-quiet or compact radio quasars, that may only have 
been active quasars for a few x 10 2 y. 



2.3 Radio size vs redshift correlation 



It has been known since the nineteen sixties that th ere is a stati stical decrease 
of the angular size of radio sources with redshift (|Milevlll968D . More recent 
work includes measure ment of the angular size-redshift relation for luminous 
exten ded radio sources (iNilsson et allll993l : lNeeser et alll995l : lDalv and Guerral 
2002) , quasars (jBuchalter et a Jl99i and compact radio sources (iGurvits et all 



1999) 



Over the years there have been many valiant attempts to use the angular- 
redshift relation to derive information about the geometry of the Universe 
(e.g. determination of go and ^o) and even to set constraints on dark energy 
(|Podariu et a]||2003h . However, there are many observational selection effects 



involved. Furthermore, it is difficult to disentangle varying geometry of the 
Universe from effects due to physical evolution of the radio sources, their host 
galaxies and the surrounding ambient medium. For example, the sizes of radio 
sources can be expected to decrease at larger redshifts due to a systematic 
increase in density of the ambient medium. Also the energy density of the 
cosmic microwave background increases as (1 + z) 4 , substantially enhancing 
inverse Compton losses. This will tend to extinguish radio sources earlier in 
their lives, making them on average smaller. 

In summary, the many non-cosmological effects that influence the angular 
size-redshift relation have made it impossible to draw robust conclusions 
about cosmology from such considerations. 



2.4 Radio spectral index vs redshift correlation 

One of the most intriguing properties of the relativistic plasma in HzRGs 
is the strong correlation that exists between the steepn ess of radio source 
spectra and the redshift of t he associated radio galaxies ( Tielens et all 1 19791: 



iBlumenthal and Milevlll979T ). Radio sources with very steep spectral indices 
at low frequencies <; 1 GHz tend to be associated with galaxies at high 
redshift (e.g. Figure [3]). This empirical correlation between radio spectral 
steepness and redshift has proved to be an efficient method for finding distant 
radio galaxies. Most known HzRGs (Table 3) have been discovered through 
following up those radio sources with the steepest ten percentile of radio 
spectra (spectral index, a ~ -1). 

The conventional explanation for the z ~ a correlation is that it is the 
result of a concave radio spectrum (see Cygnus A in Figure[3|), coupled with a 
radio K-correction. For higher redshifts, observations of sources at a fixed ob- 
served frequency will sample emission at higher rest-frame frequencies where 
the concave spectrum becomes steeper. The most important mechanisms for 
making the radio spectra concave are synchrotron and inverse Compton losses 
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Fig. 3 Left. Plot of radio spectral in dex versus redshi f t, show ing that more distant 
sources have steeper spectra. [From |Pe Br cuck ct al (|2000f )] ■ Above right. Radio 
spectrum of the luminous radio galaxy Cygnus A at z = 0.05, showing spectral 
curvature. Bottom right. Radio spectrum of the HzRG 4C41.17 and its various 
components. [From Chambe rs et all (|19901 )]. This is still one of the most well de- 
termined spectra of a HzRG at low frequencies. Note the absence of significant 
spectral curvature. 



at high frequencies (e.g. iKlamer et alH200"r3) and synchrotron self-absorption 
at low frequencies. 

Although there is evidence from radio colour-colo ur plots that systemati c 
concave curvature is present in the radio spectra (e.g. lBornancini et alll2007l ). 
such an explanation alone is insufficient to explain the observed correlation. 
The radio spectra of many distant luminous radio galaxies are not concave 
at the relevant frequencies. The radio source with the most accurately deter- 
mined spectrum over a wide frequency range is 4C 41.17 at z = 3.8 (bottom 
right in Figure [3]). This source has an extremely straight spectrum between 
40 MHz and 5 GHz, t he relevant frequency range for the z ~ a correlation 
( Chambers et all lT990). Although the spectrum steepens above 5 GHz, this 
is too high a freq uency to contr i bute to the z ~ a correlation. Furthermore, 
a recent study bv lKlamer et all ( 2006t ) showed that 33 of 37 sources in their 
SUMSS-NVSS sample have straight and not concave spectra between 0.8 and 
18 GHz. 

Two alternative effects have been proposed to explain the observed z ~ a 
relation. The first possibility is that the z ~ a relation is an indirect manifes- 
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tatio n of a luminosity, L ~ a effect ( Chambers et al 1990; Blundell and Rawlingsl 
Il999h . 

Classical synchrotron theory predicts that continuous particle injection 
will result in a spectrum with a low- frequency cut-off, vi, whose frequency 
depends on the source luminosity, L according to v\ ~ L~ 6 / 7 . For a flux- 
limited sample, Malmquist bias will cause sources at higher redshift to have 
preferentially larger radio luminosities. Over the relevant frequency range, 
the L -6 / 7 vs vi effect would therefore result in the observation of an z ~ 
a relation. However, it is unlikely that this luminosit y - spectrum relation 
i s the correct explanation, or at least the whole story. lAthreva and Kapahil 
(1998) showed that a z ~ a correlation still persists even when samples are 
restricted to a limited range of L. 

A second explanation is that some physical effect causes the spectral 
index to steepen with h igher ambient density and that the ambient de nsity 
increases with redshift (jAthreva and Kapahil [19981 iKlamer et alll2006l ). For 
example, putting the radio source in a denser environment would cause the 
upstream fluid velocity of the relativistic particles to decrease and a first- 
order Fermi acceler ation process would then produce a steeper synchrotron 
spectrum. Recently IKlamer et all ( 20061 ) pointed out that such a mechanism 
would (i) result in both z ~ a and L ~ a correlations and (ii) provide a 
natural physical link between high-redshift radio galaxies and nearby cluster 
halos. 

However, it is difficult to produce the observed z ~ a relation from such a 
simple density-dependent effect alone. The clumpy UV/optical morphologies 
(e.g Section I6.1|) indicate that the density of gas around HzRGs is highly 
non-uniform and that the density is larger close to the nucleus than in the 
outer regions. Furthermore, the internal spatial variations of spectral index 
within individual source is observed to be smaller than the source t o source 
variations of the integrated spectral indices (e.g ICarilli et allll994l ). If the 
ambient medium is highly non-uniform, how can one side of a radio source 
"know" that the other side has an uncommon ultra-steep spectrum? It is 
therefore likelier that the ultra-steep spectra are produced by some mech- 
anism in which the spectral index is determined within the galaxy nucleus 
rather than by the environment at the locations of the radio lobes. 

In summary, the origin of the z ~ a effect is still unclear and more de- 
tailed information is needed about the dependence of the radio spectrum on 
redshift. In the near future accurate measurements of low- frequency radio 
spectra of HzRGs with new facilities such as LOFAR will be important for 
such studies. 



2.5 Nonthermal X-ray emission 

Although X-ray measurements of HzRGs are sensitivity-limited, significant 
progress in the field has been made in the last decade using the Chandra and 
XMM X-ray telescopes. Extended X-ray emission has been de tected from 
abou t a dozen high - redsh if t radio galax i es and radio- loud quasars (ICarilli et all 



2002t iFabian et all [200311: IScharf et all [20031: lYuan et all 120031: iBelsole et all 
2004HOverzier et al 2005; B lundell et alll2006l : lErlund et a]l200a I Johnson et all 
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2007). The extended X-ray emission is typically elongated in the direction 
of the radio source, indicating that there is some physical link between the 
X-ray emission and the relativistic plasma. 

Several mechanisms have been proposed for producing the extended X-ray 
emission. The one that is most widely sugge sted inverse Compton scatter- 
ing (IC) of the cosmic microwave background dFabian et all 20031 : IScharf et~al 
120031: IBelsole et alll2004t lErlund et all2005 1 Johnson et a]|l2007h . Because the 
density of CMB photons increase as (1 + z) 4 , IC scattering of the CMB be- 
comes increasingly important at high redshift. Under the assumption that 
the X-ray em ission is due to this pro cess, comparison of the radio and X-ray 
luminosities (jFelten and Rees) 19691) yields magnetic field strengths consis- 
tent w ith eauinartition (|Belsole et alll2004l : lOverzier et alll200"5l : I Johnson "et al 
120071 ). Because the radiative lifetimes of radio synchrotron-emitting electrons 
are shorter than the lifetimes of the X-ray emitting IC electrons, the IC emis- 
sion traces older particles. In 4 C 23.56 at z = 2.48 , the X-ray emission is 
observed to extend by ~ 500kpc ( Johnson et a]||2007t ). implying an energy in 
both relativistic and IC-emitting electrons of > 10 59 erg, an energy reservoir 
equivalent to ~ 10 8 supernovae. 

Other processes that have been invoked to produce the extended X-ray 
emission include inverse- Compton up-scattering of synchrotron photo ns in 
the jet (synchrotron self-Comp ton emission - S SC) fccharf et albooah and 
thermal emission from shocks (Carilli et a i l2002t IBelsole et alll2004D (see Sec- 
tion O)- 



3 GAS AND DUST 

"Sometimes you can't stick your head in the engine, so you have to examine 
the exhaust" Donald E. Osterbrock 

Gas and dust, in various states, are important constituents of HzRGs and 
must play a fundamental role in their evolution. Gas has been observed over 
a wide range of temperatures and having a variety of forms. 



3.1 Hot ionised gas - Radio depolarization 

X-ray observations show that hot ionized gas is widespread at the centres of 
clusters and around radio galaxies at low redshifts. Typical derived param- 
eters are temperatures of ~10 K, densities of 0.05 cm~ 2 , and masses of a 
few x 10 12 cm— 3, with pressures of ~ 10~ 9 dyne cm -2 . Such a gas would be 
sufficient to confine both the radio-emitting plasma and the Lya halos (see 
Section |3~2"|) of HzRGs. For low and intermediate redshift clusters, bubble-like 
structures in X-rays have been attribut ed to energy injection into such a hot 
intracluster medium by radio jets (e.g. iMcNamara et alll2000l ). 

The sensitivity of the present generation of X-ray telescopes is only marginally 
sufficient for detecting similar hot gas at z ~ 2. X-ray em ission has been ob - 
served from the Spiderweb Galaxy at z = 2.2 (Figured]) (jCarilli et alll2002l ). 
The X-ray emission is extended along the radio source axis and attributed to 
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Fig. 4 X-ray emission from 
the Spiderweb Galaxy, PKS 
1138-262 at z = 2.2, ob- 
served with the Ch andra X-ray 
Telesc ope [From ICarilli et all 
(2002)]. X-ray contours are su- 
perimposed on a VLA gray 
scale representation of the 5 
GHz radio continuum emission 
at 0.5" resolution. The cross 
marks the position of the radio 
galaxy nucleus. Note that the 
X-ray and radio emission are 
aligned wit h e ach other. See 
also Figure [16] 




RIGHT ASCENSION (J2000) 



thermal emission produced by shocks as the synchrotron-emitting relativistic 
jet that propagates outwards. 

An alternative technique for detecting and studying this hot gas is throug h 
its effect on the polarisation of background radio emission (e.g. lMilevl ll980). 
As it propagates through a magnetoionic medium the polarisation of linearly 
polarised radiation is rotated through an angle \ that is proportional to the 
square of the wavelength, A (cm). 

A\ = 5.73 x 10~ 3 RM A 2 , where R is the observed rotation measure, 
given by RM is RM = 812/ U t (s)S|| (s)ds - 8.1 xlO 8 n t B\\s rad m~ 2 . 

Here s(kpc) is the path length through the medium, n t (cm -3 ) is the 
density of thermal electrons and Bii is the component of magnetic field along 
the line of sight. The intrinsic rotation measure is RM = RM j, s x (1 + z) 2 

O bservations of radio polarisation data for > 40 HzRGs (ICarilli et all 

11994 Il997t lAthreva et all 1 19981 : |Pen"tericci et a]|[200fi iBroderick et all 120071 1 
revealed that several sources have extremely l arge rotation measures (RMs) , 
up to more than 1000 rad m~ 2 in the rest frame( Carilli et all 19971 : lAthreva et all 
1998) (e.g. Figure [5]), with considerable variation across the sources. An ob- 
vious interpretation of the large observed rotation measures is that HzRGs 
are embedded in dense gas. At low-redshifts the largest rotation measures are 
observed for radio sources that are located in X-ray emitting clusters that are 
inferred to be " cooling flows" (|Ge and Owenlll994l ) . Hence the large rotation 
measures of HzRGs are strong circumstantial evidence that these objects are 
also located in cluster environments at high redshifts (Section El- 
Assuming that (i) the rotation is produced over a line-of-sight path length 
that is comparable with the projected dimensions of the radio sources and (ii) 
the magnetic field strength is close to the equipartition values, the derived 
electron densities are ~ 0.05 cm -3 , comparable to those obtained for low- 
redshift clusters from the X-rays. 



3.2 Warm ionized gas 



The bright emission lines emitted by luminous radio galaxies allow their 
redshifts to be determined and have been an important reason why radio 
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Fig. 5 The observed polarisa- 
tion position angle as a func- 
tion of the square of the 
wavelength 2 for a radio knot 
of 4C41.17 at z = 3.8. [From 
ICarilli et all lfl994l )]. The knot 
is the brightest of the 4 radio 
knots shown in Fig. [5] Note 
the excellent fit to a A 2 depen- 
dance and the large resultant 
rotation measure for the knot. 



galaxies have been such important cosmological probes. The emission lines 
are emitted by warm (10 4 5 K) gas and provide a powerful tool for studying 
physical conditions within the HzRGs. 

Emission lines observed from active galactic nuclei are both "permitted" 
recombination lines and collisionally-excited "forbidden" lines. When a free 
electron is captured by an atomic nucleus it cascades to the ground state 
emitting a scries of recombination lines. The most prominent recombination 
lines are due to the most common elements, hydrogen and helium. However, 
some of the strongest lines in the spectra of emission nebulae correspond to 
so-called " forbidden" transitions that have miniscule probability of occurring 
relative to the permitted transitions, but whose energies lie within a few kT 
of the ground levels and can therefore be easily populated by collisions. Each 
element has a critical density below which most de-excitations are radiative 
and a line is produced. Forbidden lines are only produced by gas with den- 
sities in the range ~ 10 - 10 5 cm~ 3 and the line ratios contain important 
information abou t the density, temperature, ionisation and a bundances of 
the emitting gas (jOsterbrock and Ferlandl2006t |Petersodll"997| ). 

A composite emission-lin e spectrum for H zRGs, and a list of emission 
line strengths were given bv iMcCarthvl (|l993l) . During the early nineties it 
was found that, for high-redshift radio galaxies, nuclear emission lines are 
generally accompanied by an additional component that is highly extended 
spatially. The presence of giant luminous i onized gas nebulae are amongst 
the most remarkable features of HzRGs (e.g. lReuland et all2003l ). These "ha- 
los" have sizes of up to ~ 200 kpc (e.g Fig \W\i and their study provides a 
wealth of information about kinematics and physical conditions within and 
surrounding the HzRGs and on the origin of the gas. 

During the last 15 years considerable progress has been made in studying 
the properties of Lya halos, using large ground-based optical/IR telescopes. 
Observations of their emission line spectra have been extended into the in- 
frared (optical rest-frame), where additional useful physical diagnostics can 
be obtained. Also, the kinematics of several HzRGs have been mapped in 
detail. An excellent recent review of the emission-line properties of HzRGs 
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has been given bv lVillar-Martir] lj2007h . For detailed information the reader 



is referred to that review and the references therein. 

The ionised gas halos have Lya luminosities of typically ~ 10 43 5 erg s" 1 
and line widths varying from a few hundred km/s in their outer parts to 
> 1000 km s -1 , near the gala xy nuclei. The usual emission-line diagnostics 
( Osterbrock and Ferlandll2006l ) show that the gas has temperatures of T e ~ 



10 4 K - 10 5 K, densities of n e - 10 5 - 10 15 cm" 3 and masses of ~ 10 9 - 
10 10 M Q . The warm gas occupies a relatively small fraction of the total vol- 
ume of the HzRGs, with filling factors estimated to be ~ 10~ 5 compared 
with unity for the hot gas. Although the topology of the gas is not well 
determined, the properties of the emission lines and the covering factor de- 
duced from statistics of absorption lines (Section l3.3.2p clouds led to a model 
van Oiik et allfl997l l in which the halos are composed of ~ 10 12 clouds, each 



having a size of about 40 lig ht days, i.e. comparable with that of the solar 
system. Ivan Oiik et all ( 19971 ) speculate that the clouds might be associated 



with the early formation stages of individual stars. 

The morphologies of the Lya halos are clumpy and irregular. Their overall 
structures are often aligned with the radio axes and sometimes extend beyond 
the extremities of the radio source. There appear to be two distinct regimes 
in the halos that often blend together. The inner regions close to the radio 
jets are clumpy, with velocity spreads of > 1000 km s . They appear to 
have been perturbed by the jets. The outer regions are more quiescent, with 
velocity ranges of a few hundred km s^ 1 ). They appear to be more relaxed 
than the inner regions (e.g. Fig. [TJ). 

The relative intensities of emission lines are, in principle, powerful diag- 
nostic tools for studying physical conditions in the warm line-emitting gas. 
However, disentangling the effects of ionisation, abundances, density and tem- 
perature using the emission line ratios is complicated and requires detailed 
modeling, incorporating all facets of the HzRGs. We shall now discuss the 
physical conditions within the warm gas halos in more detail. 

3.2.1 Ionisation 

Various mechanisms have been proposed for exciting the gas. These include 
(i) photoionisation from an AGN, (ii) photoionisation from stars, (iii) pho- 
toionisation by ionizing X-rays emitted by shocked hot gas and (iv) col- 
lisional ionization from shocks. Plots of optical- line ratios have been used 
extensively to study the ionisation of gas in nearby active galaxies, where 
evidence for both jet- and accretion-powere d shocks and for photoionization 
by the central AG N has been found (e.g. IVillar-Martin and Binettd Il997t 
Villar -Martin et a 1 119971 : iBicknell et aS oOO: I Groves et a 1l l2004l f). 



Although these relationships have been most accurately calibrated in the 
rest-frame optical region of the spectru m, line diagnostic diagrams have als o 
been developed for use in the ultraviolet (|Allen et allll998tlGroves et alll2004 . 



the spectral region of HzRGs sampled by optical observations. Interpreta- 
tion of emission line ratios are complic ated by effects of dust and viewing 
angles (IVillar -Martin and Binettelll996| ). Also most HzRG spectra are spa- 
tially integrated over regions of 10 or more kiloparsec, where conditions and 
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sources of excitation may change. In a recent comprehensive study using as 
ma ny as 35 emissio n lines throughout the rest-frame UV and optical spec- 
tra Humphrey et al (2006) concluded that photoionization is the dominant 
source of excitation in the quiescent gas. A harder source of photoionisation 
than stars is needed, consistent with photons from an AGN. 

If a quasar is exciting the warm gas in the halo, why don't we see it? 
The usual explanation is that the quasar emits radiation anisotropically. It 
is highly absorbed in the direction in which we view it, but not along the 
radio axis. Support for this idea is provided by the large optical polarisations 
measured for some HzRGs. (Section 13.5. 2|) . However, if highly anisotropic 
flux from a hidden quasar is indeed the dominant source of excitation, it 
is difficult to understand why many of the Lya halos are approximately 
symmetric in shape. An alternative explanation is that the quasar activity is 
isotropic and highly variable, with short sharp periods of intense activity and 
longer periods of relative passivity, when the quasar is dormant (see Section 

E3J. 

Although photoionisation by a quasar is presently the " best bet" for the 
dominant source of excitation, it is not likely to be the only culprit. Con- 
siderable variation is observed in the relative emission line strengths from 
object to object and within individual HzRGs. There is strong evidence that 
there is also c qllisional excitation from shocks, par ticularly close to the radio 
jets - see (e.g. iBicknell et a il l200nt iBest et all 12000) and Section l4~3l 



3.2.2 Abundances and star formation 



Although Lya is by far the brightest emission line emitted by the warm 
gas halos, other lines are also detected. In 4C41.17 at z = 3.8, kinematical 
structures in the Lya line are closely followed in the carbon CIV and oxygen 
[Oil] and [OIII] lines, with the [QUI] A5007 em ission extending by as much 

as 60 kpc from the nucleus (jReuland et a]|[2007l ). See also Fig. [7] 

T he chemical abundance of the halo gas is close to solar ( Vernet et all 

120011: iHumphrev etai [20061. consistent with the HzRG having undergone 
prodigious star formation at earlier epochs. Further indications that the star 
formation rate in HzRGs was higher in the past come from measurements of 
the Ly a luminosities and Lya/Hell ratios that are both s ystematically larger 
for H zRGs with z > 3, than for those with 2 < z < 3 ( Villar-Martm et all 

Hoi). 

The relative intensity of the NV A1240A line varies from being an order o f 
magnitude fainter than the carbon and helium lines dDe" Breuck et a 1 l2000h 
to being as luminous as Lya ( van Oiik et al|[l994L e.g.). This has been in- 
terpreted as indicating that there are large variations in metallicity between 
HzRGs. 



3.2.3 Inner Kinematics: Outflows, jet interactions and superwinds 

In considering the kinematics of the ionized gas halos, we shall separately 
consider the nuclear turbulent regions, where there is evidence for outflows 
and the outer passive regions where the most likely dominant systematic 
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motions are inwards. The inner regions show shocked gas closely associated 
with the radio lobes. These display disturbed kinematics and have expansion 
velocities and/or velocity dispersions > 1000 km s _1 . 

Besides synchrotron jets there i s evidence that starburst "superwinds" 
( Armus et all 1 19901 IZirm et al l2005h are also present in the inner regions of 



the Lya halos. Both the jets and the superwinds will exert sufficient pres- 
sure on the warm gas to drive it outwards. Interactions between radio jets 
and the ambient gas is important in some low-redshift radio galaxies (e.g. 
iHeckman et~alll982l Il985 Ivan Breugel et allll984l Il985l |). The radio sources 
are observed to excite, disturb and entrain the gas. Likewise, the gas can 
bend and decollimate radio jets and enhance the intensity of their radio 
emission through shock-driven particle acceleration. In general small r adio 
sources show more jet-gas interactions than large ones (|Best et alll2000h . At 
z > 2, HzRG exhibit signatures of even more vigorous jet-gas interactions 
i.g. Villar-Martin et a 1 1998; Villar -Martin et a ll999Ll200afHumphrev et~all 



2006). The kinematics is more turbulent and the ionisation is higher in the 
reg ion of the jet than in the quiescent outer halo ([Humphrey et alll2006l ). 
iNesvadba et all ( 2006h recently carried out an important study of the Spi- 



derweb Galaxy with an integral field spectrograph in the optical rest frame. 
They make a convincing case that there is an accelerated outflow of warm gas 
in this object. The only plausible source of energy for powering this outflow 
(few x 10 60 erg) is the radio jet and, even then, the coupling between the jet 
and the ISM must be very efficient to account for the observed kinematics. 
The pressure in the radio jets can drive gas outwards from the nuclei for tens 
of kiloparsecs and play an important role, together with starburst-driven 
superwinds, in "polluting" the intergalactic medium with metals. There is 
spectroscopic evidence for the ejection of enriched mat erial in 4C 41.17 at z 
= 3.8 up to a distance of 60 kpc along the radio axis (jReuland et all 12007t h 
See also Fig. [7] 

In Section l4~3l we shall discuss strong evidence that the jet-gas interactions 
can also trigger star formation in HzRGs. 

3.2-4 Outer kinematics: Infall of the quiescent halos 

An important diagnostic in tracing the origin of the warm gas is the kine- 
matics in the outer region of the giant halos, where they are apparently un- 
perturbed by the radio jet. The outer halos displays systematic velocity vari- 
ations of a few hundred km/s (e.g. s ee Figure E). Are these systematic veloc- 



ity variatio ns the resu lt of rotation (Ivan Oiik et all 19961: IvTllar-Martm et all 
12001 120061) . outflows (|Zirm et a 1 2005; INesvadba et alll2006| ) or infalling mo- 



tions? 

It is difficult to discriminate between the various kinematic scenarios from 
velocity data alone, but comparison of spectroscopic and radio data pro- 
vides additional relevant information. In a study of 11 HzRGs with redshifts 
2.3 < z < 3.6, several correlated asymmetries were f ound between the halo 
kinematics and assymet ries in the radio structures ( Humphrey et all 120071 : 



Villar -Martin et a fl l2007f) . On the side of the brightest radio jet and hot spot, 
the quiescent nebula appears systematically redshiftcd (receding) compared 
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Fig. 6 Velocities and velocity dispersi ons across MRC 2104 -242 at z = 2.5 mea- 
sured with VIMOS on the VLT. [From IVillar-Martfnl (120071 )]. This HzRG is sur- 
rounded by a giant Lya nebula that extends by ~ 120 kpc The position of the radio 
core is indicated with a cross. The velocity field a ppears symmetric an d ordered 
implying either rotation or radial motions. See also Ivan Oiik etlaj {1996). 



with the other side. On the assumption that the bright radio j et and hot 
spot are moving t owards us and brightened by Doppler boosting ( Reesl[l"967t 
lKellermannll2003| ) . the quiescent gas must be moving inwards. The brightest 
radio hotspot is also the least depolarised one, as expected if it is on the 
closest side of the HzRG and consistent with the infall scenario. 

Could the inflowing motion of the gas be a result of cooling flo ws, that 
have long been studied in c lusters of galaxies at low redshifts, (e.g. IFabianI 
Il994t iKaufmann et all l20061? Recent XMM and Chandra observations have 
shown that the cooling rates are reduced by an order of magnitude below the 
simple cooling flow models at temperatures < xl0 7 K, probably through in- 
teraction of the gas with radio sources associated with galaxies in the cluster. 
An alternative diagnostic of cooling flows that are more feasible for HzRGs 
than X-rays is the measurement of roto-vibrational lines from H2 molecules 
at ~ 2000K (jJaffe and Bremerlll997ilJaffe et alll2005l ). At high redshifts these 
lines are shifted from the near to mid-IR bands, accessible with the Spitzer 
Telescope. The H2 lines are a promising tool for studying both the cooling 
in the gas around HzRGs and the excitation mechanisms. 

To summarise, the kinematics of the warm gas is complex. There is evi- 
dence that gas in the outer regions of the halo is flowing inwards, providing 
a source of material for feeding the active nucleus (Section [5]) and that gas 
in the inner region is being driven outwards by pressure from jets and star- 
bursts. The various motions in the ionized gas halos are likely to contribute 
to feedback processes between the AGNs and the galaxies invoked in current 
models of massive galaxy evolution (Section [5]) . 
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Fig. 7 L ong slit spectroscopy of the ionized gas halo in MRC 0943-242 at z = 
2.9. [From Villar-Marti'n (|2007j ): iHumphrev et all I|2006I 1]. VLA radio map overlaid 
with the WFPC2 HST image (left), spatially aligned with the 2-dimensional Keck 
spectra of the main UV rest frame emission lines. The impact of jetgas interactions 
on the observed properties of the giant nebula (Size ~ 70 kpc) can be seen in the 
much broader and brighter emission lines within the radio structures compared 
with the faint emission detected beyond. 



We end our discussion of the halo kinematics on a cautionary note. Most 
kinematic studies of HzRGs until now have been based on Lya, because of 
its relatively large equivalent width and its accessibility with optical tele- 
scopes. However, Lya is a resonant line and subject to strong scattering and 
optical depth effects. Hence, the resultant kinematics may not be completely 
representative of the g whole, particularly in the inner regions. 

3.2.5 Relation to non-radio Lya nebulae 



There may well be a connection between Lya halos in HzRGs and the dis- 
emb odied extended Lya nebulae that have been discovered i n recent years 



(e.glFvnbo et allll99a ISteidel et all I2Q00L iFrancis et alll200ll : iMatsuda et~al 
2004 iDev et all 120051 : IColbert et all 120061 : iNilsson et alll2006h These nebulae 



[2004 

also have physical extents ;> 100 kpc and Lya line fluxes of ~ 10~ 15 ergs s _1 
cm -2 . Although in many respects, these Lya blobs resemble the giant ionised 
gas halos around HzRGs, they have < 1% of the associated radio continuum 
flux and no obvious source of UV photons bright enough to excite the nebu- 
lar emission. However, millimeter emission has been detected from seve ral of 
these nebulae dChapman et alll200ll : ISmail et alll2003l : iGeach et alll2005l ) and 
IMatsuda et all ( 2004 ) suggest that the extended Lya nebulae are also associ- 
ated with dense environments in the early Universe. Just like the HzRG halos, 
these nebula e can be excited by qua s ars that are heavily obscur ed along the 
line of sig ht (|Haiman and Reesil200ll : IWeidinger et~all2004 l2005t l , or quasars 
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that undergo recurrant flares (Section 15. Alternatively, they may be asso- 
ciate d with cooling-flow-like phenomena ( Haiman et all BoOOt IDijkstra et all 

Hoi). 



3.2.6 Nebular continuum 



iDickson et all ( 19951 ) pointed out that the nebular continuum emission due 



to the line-emitting warm gas of HzRGs is significant and must be taken into 
account when computing the various contributions to the U V continuum. 
See also (|Allerlll984t IVernet et alll200li iHumphrev et~all2007l ). However, in 



general the nebular component contributes < 25% of the continuum emission 
at 1500 A, and is much less important than starlight or emission from a 
hidden or dormant quasar. The contibution from nebular continuum can be 
quite accurately predicted by means of the st rength of t he Hell A4686A line, 
or when not available the Hell A 1640 A line (|Allerlll984l ). 



3.3 Neutral Gas 

There are two techniques for studying cool HI gas in HzRGs. One method 
is to measure redshifted absorption of the 21 cm hydrogen line against the 
bright radio continuum. The second method is to observe deep narrow ab- 
sorption troughs that are often present in the Lya profiles. In principle, both 
these techniques can be used to constrain properties of the neutral hydrogen 
such as spatial scales, mass, filling factor, spin temperature and kinematics 
(|Rottgering et allll999h . 



3.3.1 HI Absorption 

Neutral hydrogen (HI) atoms are abundant and ubiquitous in low-density re- 
gions of the ISM. They are detectable by means of the hyperfine transition, 
emitting at 1420.405751 MHz (~ 21 cm). Observations of atomic neutral 
hydrogen by means of this line has been one of the most powerful tools of 
radio astronomy since its inception. Studying the 21cm HI line in absorp- 
tion has been an import ant probe of HI around low-redshift radio galaxies 
see the recent review of iMoreantil (|2006h . Limitations of such searches for 



HI in high-redshift radio galaxies include the availability of low-noise re- 
ceivers that cover the observing frequencies dictated by the target redshifts 
and problems due to radio frequency interference. In 1991 HI absor ption 
was detected in the HzRG, 0902+34 at z = 3.4, bv lUson et al l| 1991ft . See 



Figure [51 Since then prog r ess in this field has been disappointing. For a re- 
view see iRottgering et all (1 1999ft. Besides some follow-up work on 0902+34 



view see nottgcrmg ct al (lyyy). rjcside s some rollow-up work on uyu^+d4 
(|Briggs et aJl993tlDe Bruvnll996HCodv and Braud200alChandra et alll2004h . 



there have been only two tentative but unconfirmed additional detections of 
HI absorption in 0731+438 at z = 2.4 and 1019+053 at z = 2.8. 

The HI absorption line provides a measure of the average column density 
of the absorbing material weighted by the flux density of the background 
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Fig. 8 Absorbing neutral 
hydrogen. Top: The red- 
shifted HI absorption profile 
of 0902+34 at z = 2.4. [Fro m 
ICodv and Braunl ([20031 )]. 
Bottom: The redshifted Lya 
absorption profile of 0943242 
at z = 2.9, wit h a model 
overla id. [From Jar vis et all 
(2003)]. 
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source. The total colum n density i s ~ 4. 4 x 10 18 T5. Assuming a spin tem- 
perature of T s ~ 10 3 K. lDe Bruvr] (|l996t ) derives an HI column density of ~ 



2 x 10 21 cm 2 for 0902+34 and a corresponding mass of neutral hydrogen 
ofM HI ~ 3 x 10 7 M Q . 

Why has HI absorption not been detected in a larger number of HzRGs, 
despite extensive searches? A possible explanation is that the HI absorption 
is caused by small ~ 100 pc-si zed disks or torus-like structures, aligned per- 
pendicular to the radio source! Rottgering et all fl999). The radio emission of 
0902+34 is more centrally concentrated than that for most HzRGs, consis- 
tent with the hypothesis that the absorption is produced by such a small 
disk. For the more extended radio sources associated with most HzRGs, the 
disk covering fraction would be very small and the disk would not produce 
significant absorption. 



3.3.2 Lya absorption. 



More than a decade ago lvan Oiik et al ( 19971 ) discovered that strong absorp- 
tion features are common in the Lya profiles of HzRGs. Such features were 
present in the majority of the 18 HzRGs that they studied with sufficient 
spectral resolution. Derived column densities were in the range 10 18 — 10 19 5 



cm 



The absorption, usually interpreted as being due to HI surrounding the 
HzRG, provides an interesting diagnostic tool for studying and spatially re- 
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solving neutral gas surrounding HzRGs. Because the spatial extension of 
the absorbing region can be constrained, the Lya absorption lines provide 
information about properties of the absorbing gas (e.g. dynamics and mor- 
phologies) that cannot be studied using quasar absorption lines. Since in 
most cases the Lya emission is absorbed over the entire spatial extent (up 
to 50 kpc), the absorbers must have a covering fraction close to unity. From 
the column densities and spatial scales of the absorbing clouds, the derived 
masses of neutral hydrogen are typically ~ 10 s M@. 

Additi onal information abo ut the properties of the HI absorbers was ob- 
tained bv IWilman et all (|2004l ). In a study of 7 HzRGs with 2.5 < z < 4.1, 
they identified two distinct groups of H I absorbers: strong absorbers with 
column densities of N#j ~ 10 18 to 10 20 cm -2 and weaker systems with N#j 
~ 10 13 10 15 cm~ 2 . They suggest that the strong absorbers may be due to 
material cooling behind the expanding bow shock of the radio jet and that 
the weak absorbers form part of the multi phase proto-in tracluster medium 
responsible for the Lya forest. Furhermore. iKrausd ( 2005h carried out hydro- 
dynamic simulations of a HzRG jet inside a galactic wind shell and showed 
that strong HI absorption could be produced. 

Absorption is occasionally observed in the pr o files of other emiss ion lines 
than Lya, such as CIV - ([Rottgering et aliri995l ). iJarvis et~al (|2003h studied 
the profiles of two of the most prominent absorbing HzRGs 0943-242 at z 
= 2.9 and 0200+015 at z = 2.2 with high spectral resolution. The data are 
consistent with a picture in which the absorbing gas has low density and low 
metallicity and is distributed in a smooth absorbing shell located beyond the 
emission-line gas. However, the metallicity, inferred from the C IV absorption, 
is considerably lower in 0943-242 than in the slightly larger source 0200+015. 
This difference in metallicity is explained as due to chemical enrichment 
via a starburst-driven superwind (Section 13.2.311. Further o bservations and 
modeling of the spectrum of 0943-242 bv lBinette et all ( 2000h in dicate that in 
this o bject the absorbing gas may actually be ionised - see also (jBinette et all 
120061 ). However, 0943-242 may be a special case. It has a relatively small radio 
size and one of the deepest Lya absorption trough of all known absorbers. 



3.4 Molecular Gas 

Star formation is generally observed to occur in molecular clouds - cold dark 
condensations of molecular gas and dust, that are observable in the millime- 
ter and near-IR. In these clouds atomic hydrogen associates into molecular 
hydrogen, H2, a species that unfortunately does not emit easily observable 
spectral lines. The next most abundant molecule is carbon monoxide. Because 
the rotational transitions of the dipolar molecule 12 CO are caused primar- 
ily by collisions with H2, CO is an excellent tracer of molecular hydrogen. 
The most important redshifted CO transitions for the study of high-redshift 
objects are J=(l-0), (2-1), (3-2), (4-3) and (5-4) at 115.2712, 230.5380, 
345.7960, 461.0407 and 576.2677GHz respectively. These lines are an im- 
portant diagnostic for probing the reservoir of cold gas available for star 
formation. 
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Table 2: HZRGs detected in CO 
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1.2±0.3 




9 



"Only broad component; narrow component has 5*00^^=0.09+0.01 Jy km s 
Refer ences: 1 =IScoville et all (IT997T). 2 = IAlloin et all (|2000T), 3 = IDe Breuck et all 
(l200l). 4 = IDe Breuck et all (120031). 5 = IPapadopoulos et all (12000 ) 6 = 
IGreve et all (120041 ) , 7 = IPapadopoulos et all (|2005l ). 8 = IDe Breuck et all ((2003) . 
9 = lKlainCTet "aT (l2005l ). 



Intensive searc hes for CO emission from HzRGs du ring the early 1990s 
were unsuccessful (|Evans et a 1 199(5; van Oiik et alll997l) . Since then, the sen- 
sitivity of (sub)millimctrc receivers has been improved and the high-redshift 
Universe has been opened to molecular line studies (see reviews by Solomon 
& vanden Bout (2005) and Omont (2007)). 

To convert the CO to H2 mass, one often assumes a standard conversion 
factor. For high redshift CO studies, this factor is calibrated b ased on obser- 
vatio ns of nearby ultra- luminous infrared galaxies (ULIRGS) (jDownes et all 
Il993| ). With the assumption that this value is also applicable to high redshift 
objects, one can use the strength of the (1-0) CO transition to derive the 
mass of the molecular gas. Because current centimetre telescopes do not allow 
observations of the (1-0) transition at z < 3.6, one needs to observe higher 
order transitions shifted to the atmospheric windows at 3, 2 and 1.3 mm. 

The inferred masses of H2 in the CO-detected galaxies are between 10 10 
and 10 11 Mq, indicating that there is a large mass of molecular gas in these 
objects and a substantial reservoir of material available for future star forma- 
tion. However, the calculated masses should be treated with ca ution, because 
their derivation is based on a large number of assumptions (e.g. lDownes et all 
Il993l ). Observations of higher CO transitions are biased to the detection of 
denser gas than ground-state observations and can result in an underestimate 
of the total molecular gas mass. 

Observations of multiple rotational transitions allow the temperature and 
density of the molecular gas to be constrained using large velocity gradient 
models. The results indicate that the CO properties are heterogeneous. TN 
J0924-2201 is only detected in the ground-state (1-0) transition, but 4C 41.17 
i s only detected in CO(4-3), despite sensitive searches for the (1-0) transition 
( Ivison et a 3 119961: IPapadopoulos et"al l2005t l. This high excitation level in 
4C41.17 implies large gas densities n(H_2) > 1000 cm -3 , consistent with gas 
fueling a nuclear starburst. 

In a few cases the CO appears spatially resolved ( Papadopoulos et "al 

12000b IGreve et all 120041; IDe Breuck et alll2005D . extending over 10 - 20 kpc 
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(e.g. see Figure El), pro viding; kinematic information about the molecular 
gas. lKlamer et alT(|2004l ) have claimed that there are alignments between the 
molecular gas and radio morphologies in some of the detected CO HzRGs, as 
might be expected from jet-induced star formation (Section I4.3|) . However, 
higher resolution, larger signal to noise and more statistics are needed before 
any conclusions about possible alignments can be drawn. 

Offsets between the velocities of the molecular gas (CO) and those of 
the warm gas (e.g. H e IIA16 40) of up to 500 km/s have been measured 
(|De Breucket all l2003l l. l2005h . because the H? masses exceed those of the 



warm ionized gas Lya halos by an order of magnitude, the CO lines provide 
a better measure of the systemic redshift of the HzRGs than the UV and 
optical emission lines. 

The CO detections listed in Table l3~4l are not representative of the distri- 
bution of CO in HzRGs. There are several observational biases inherent in 
current CO studies: 

— Many searches for CO at high redshifts pre-select targets on the basis 
of prior dust detections, i.e continuum millimetre flux (see Section [3.5. 1| . 
However, 850/im continuum emission was not detected for two of the CO- 
emitt ing HzRGs - 53W002 at z = 2.3 90 (IScoville et allll997l : lAlloin et all 



\200d } and TN J0924-2201 at z = 5.2 (jKlamer et alll2004h . This indicates 
that the gas/dust ratios in high redshift CO samples could be be system- 
atically underestimated due to bias in favour of dusty targets. 

— Positive detections are more likely to be published than non-detections. 
Furthermore, the derived limits are often dubious in view of relatively 
large uncertainties in the targeted redshifts. These uncertainties are due to 
large offsets that are often present between the CO and optical redshifts, 
the unknown widths of the CO emission lines, and the limited bandwidth 
of correlators used for the detections. 

— For higher redshift objects progressively higher CO transitions need to be 
observed consistent with the atmospheric transmission bands acessible to 
millimetre arrays. Hence, measurements of CO in will tend to trace denser 
gas in HzRGs than at lower redshifts. However for z > 3.6, the ground- 
state transition becomes observable with centimetre facilities, such as the 
VLA and ATCA, resulting in a discontinuity in the studies of multiple 
CO transitions needed to constrain the densities and temperatures of the 
gas. 

To investigate the statistics of CO in HzRGs and the evolution of molec- 
ular gas, it is important to carry out deep CO observations on well defined 
samples of HzRGs in multiple CO transitions, and without a pre-selection on 
their cold dust properties. The faintness of the CO lines and the uncertain- 
ties in their redshifts means that comprehensive searches for CO presently 
require large amounts of observing time. Limiting factors are the sensitivi- 
ties of current facilities and the relatively limited instantaneous bandwidth 
covered by the receiver back-ends. 

Opportunities for observing molecular gas in HzRGs will be revolutionised 
during the next decade by ALMA and the EVLA. These facilities will pro- 
vide an enormous improvement in sensitivity and discovery space throughout 



28 



George Miley, Carlos De Breuck 



Fig. 9 CO in 4C41.17 at z = 3.8 
with the Pl ateau de Bure int e rferom - 
eter. [From [DeBreuck et all (2005)]. 
Shown is the velocity plotted against 
position offset of CO (4-3), extracted 
along a PA of 51° (see radio mor- 
phology in Figure [2}. The central fre- 
quency at 96.093 GHz (z =3.79786) 
is based on the wavelength of the 
optical HeIIA1640Aemission line. The 
position-velocity diagram shows the 
two components of the CO emission. 
One of these is coincident with the ra- 
dio core, while a fainter component 
is spatially offset from it towards the 
southwestern lobe A. 
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the sub- millimeter, millimeter and centimetre wavebands. Furthermore, their 
resolution will allow the spatial distribution of the molecular gas in HzRGs to 
be studied. On a shorter timescale, the much wider tuning range and broad- 
band correlators that are presently coming on-line (e.g. at the EVLA and the 
ATCA) will provide new opportunities for investigating CO in high-redshift 
objects, particularly the ground state transition. 

The increased sensitivity of ALMA will also open the study of fainter 
molecular and atomic lines in HzRGs. Lines such as HCN or HCO + probe 
at least an order of magnitude denser molecular gas than CO, and are 
therefore better trac ers of the dense molecular cores in star-forming regions 
(|Papadopoulosll2007l) . 



3.5 Dust 



Dust is an important constituent of HzRGs and an additional diagnostic 
of star formation. It is both a major constituent of the molecular clouds 
from which stars generally form and an indicator that substantial star for- 
mation has already occurred. The presence of dust means that chemically 
enriched material is present. Thermal re-radiation from dust often dominates 
the spectral energy distribution of HzRGs at millimetre and sub-millimetre 
wavelengths (Fig [J) and dust is observable as a polarizing and absorbing 
medium in the optical and ultraviolet. 



3.5.1 Millimeter and sub-millimeter emission 



Since the detectio n of 4C41.17 at z = 3.8 (|Dunlop et all994lChini and Kruegell 
Il994t llvison|[T995 ) a large number of HzRGs have been o bserved at mill i meter 
wavelengths (| Archibald et afefJOli iReuland et al2004h . iReuland et al (|2004h 
analysed a sample of 69 radio galaxies with 1 < z < 5, detected at 850 and/or 
450 /im. Isothermal fits to the submillimetr e spectra give dust m asses of a 
few x 10 8 Mq at temperatures of ~ 50 deg ([Archibald et alll200ll ). 
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Possible sources of heating for the dust are X-rays from an AGN (Section 
|5|) and UV photons from young stars. The typical submillimetre luminosity 
(and hence dust mass) of HzRGs strongly increases with redshift, with a 
(1+z) 3 dependence out to z ~ 4. This is evidence that star formation rates 
were higher and/or the quasars brighter in HzRGs with z > 3. 

Disentangling whether the dust is heated by an AGN or young stars is 
difficult, because the two processes produce dust at partially overlapping 
temperatures. Evidence for heating by young stars comes from the slope of 
the Rayleigh- Jeans part of the thermal dust emission - observed at millime- 
tre and submillimetre wavelengths. This indicates that the dust tempera- 
tures iare_j^tatiyehjf cjjc^ emission 
(e.g. lArchibald et alll200lt I Stevens et alll2003tlReuland et alll2004l ). However, 
there is also evidence for a warm er(~ 300K dust) component in radio galax - 
ies, consistent with AGN heating ( Rocca-Volmerange and Remazeil lcs 2005). 



Furthermore, recent Spitzer observations at 5 < A b s < 70 fim showed that 
the Wien tail part of the dust SED is inconsistent with low dust tempera- 
tures. Dust with temperatures of ~ 300K or higher is required, consistent 
with AGN heating. Although the submillimetre luminosity of the HzRGs 
is uncorrelated with radio luminosity, sugge sting no strong depe ndence on 
the strength of the quasar/ AGN emission. (jReuland et all 12004 ) . the rest- 
frame 5 iim luminosi ty does appear to be correlated with radio luminosity 
((Seymour et alll2007h . 



It is likely that young stars and AGN both play a role in heating the 
dust. Clarification of their relative importance await accurate measurement 
of the SED near the peak of the thermal dust emission by the Herschel 
Telescope and high spatial resolution observations of the radial profiles of 
the thermal dust emission with ALMA. Once properly isolated, starburst- 
heated dust emission will be a powerful tool for measuring the star-formation 
rates in radio galaxies. On the assumption that most of the submillimetre 
emission is heated by stars, derived star formation rates are up to a few 
thousand M yr , consistent with those obtained from UV absorption line 
measurements (Section 14.21) . 

3.5.2 UV continuum polarisation 

The rest- frame UV polarisation is an important probe of dust in the inner 
regions of HzRGs (Idi Serego Alighieri etal Il989l : Idi Serego Alighieril Il997t 
ICimatti et a1ll993lll998l : IVernet et afeOoij T The fractional polarisation pro- 



vides a unique tool for determining the contribution and nature of of scat- 
tered light, while the polarisation angle allows the location of the scattering 
medium to be pinpointed. However, polarisation measurements of HzRGs 

require both high sensitivity and high precision . 

Eight of 10 HzRGs studied by lVernet et all ( 200ll ).show high continuum 



polarization just redward of Lya, with fractions, f p ranging from 6% to 20% 
(e.g. Fig. [TUJ). The shape of the polarised flux (= percentage polarisation 
x total intensity) is very similar to that of an unobscured quasar (Type 1 
AGN) in both the slope of t he continuum and the presence of broad emis- 
sion lines (|Vernet et alll200lh . This is strong evidence for the presence of an 
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Fig. 10 Spectral and polar- 
ization properties of 0211- 
122 at z = 2.3 taken from 
observations wit h the Keck 
Telesc ope. [From IVernet et all 
pOOTV. Plotted from top to 
bottom are: the observed to- 
tal flux spectrum in units of 
10" 17 erg s" 1 cm" 2 A" 1 at 
two different scales, the first to 
show strong emission lines and 
the second to show the contin- 
uum, the crosses respectively 
indicate continuum and nar- 
row emission lines (including 
their underlying continuum) . 
Shaded areas indicate regions 
of bright sky emission. 
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obscured quasar in the nuclei of HzRGs. The spectral energy distribution 
of a typical quasar is bluer than most other components that contribute to 
the SED of HzRGs between the UV and near-IR. The dilution by the other 
components is smaller in the rest-frame UV. At z > 2, this waveband is 
redshifted into the optical, where the most sensitive (spectro-)polarimeters 
exist on large telescopes, allowing detailed studies despite the cosmological 
brightness dimming. 

The polarisation angle is in most cases closely perpendicular to the radio 
structure, indicat ing that the scatt ering occurs within the cones traced out 
by the radio jets ( Vernet et a]|l200lh . The spatially resolved narrow emission 
lines (see Section [3~2)l are generally not polarised (e.g. Fig. [TO]) , indicating 
that the scattering medium must be located between the broad and narrow 
line regions. 

The wavelength dependence of the polarised emission is consistent with 
both grey dust or electron scattering. However, electron scattering can be 
excluded as it is much less efficient and would i mply masses of the scattering 
medi um which are close to the total galaxy mass ([Manzini and di Serego Alighieril 
1996). A clumpy scattering dust medium can also produce the observed con- 
tinuum slope. Using spectropolarimetric observat ions near the rest frame 
2200 A dust feature, ISolorzano-Iharrea et all (|2004j ) argue that the composi- 
tion of the dust is similar to that of Galactic dust. 
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ICimatti et all ( 1993h presented a dust-scattering model that explains the 
structure, the polarization properties and the spectral energy distribution 
of the ultraviolet aligned light with optically thin Mie scattering of quasar 
radiation emitted in a cone of ~ 90°opening angle. The required amount of 
spherically distributed dust is (1-3) x 10 8 M Q , consistent with the estimate 
for the dust mass fr om the sub- millime ter data (Section l3.5.1|) . However, more 
recent predictions ( Vernet et all l200lh indicate that the dust responsible for 
scattering the AGN emission should have a sub-millimetre emission that is 
an order of magnitude fainter than observed. 

Not all HzRGs hav e significant polarisation. Two of eight objects studied 
by IVernet et a l (2001) had fractional polarisation f p < 2.4%. There is no 
evidence that the nuclei of these HzRGs contain obscured quasars (Section 
5.ip . We discuss information derived from the UV polarisation observations 
about the obscured or dormant quasar in Section 15. II 



4 STARS 

Captain Boyle: "I ofen looked up at the sky an' assed meselj the question - 
what is the stars, what is the stars?" 

J oxer Daly: "Ah, that's the question, that's the question - what is the stars?" 
Act 1 "Juno and the Paycock", Sean O'Casey, Irish playright (1924) 

During the last fifteen years there has been mounting evidence that stars 
are a major constituent of HzRGs. Stellar signatures include (i) a plateau 
in the spectral energy distribut ion (SEP) between 1 and 2 /xm characteristic 
of sta rs - the "1.6 ^m bump" ( Simpson and Eisenhardtl[i999HSevmour et all 



l2007f ) (Fig. [T]) , (ii) the detecti on of rest-frame UV stellar absorption lines in 

a few HzRGs (|Dev et allll997l ) and (iii) measurement with the H ST that the 

morphologies of several HzRGs follow de Vau couleurs profiles, (e.g. Ivan Breugel et all 
ll998tlPentericci et alll200l IZirm et alll2005l ). 

In principle, spectral energy distributions (SEDs) are a powerful diagnos- 
tic of stellar populations and galaxy evolution. Young (<0.5 Gyr) stars domi- 
nate the SED in the UV, whereas the SED longward of 4000A gives informa- 
tion about populations older than ~ 1 Gyr. However, disentangling the stellar 
contribution to the SEDs from the nonstellar components is complicated (see 
Fig. [TJ. At rest-frame UV wavelengths, contributions from scattered quasar 
light and nebular continuum must be taken into account, whereas re-radiation 
by dust contributes to the SED at far infrared and millimeter wavebands (e.g 
Fig. 2). 



4.1 Old stars: HzRGs as the most massive galaxies 

The SEDs of old (>1 Gyr) stellar populations peak in the near-IR. The galax- 
ies with the largest K-band luminosities in the early Universe are HzRGs. 
Because old stars are the best tracers of stellar mass, the brightness of the 
observed K— band emission has long been used to ar gue that distant radio 
galaxies are very massive, with masses of up to lO 12 M0( Rocca-Vofmerange et all 
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Fig. 11 The IR SED of 

4C40.36 at z = 2.265, with 
modeling, using data taken 
with IRAC and MIPS on 
the Spitzer Te le scope . [From 
ISevmour etU lf2007t )]. The 
solid dark line indicates 
the total best-fit SED. The 
stellar and dust components 
are indicated by dashed 
and dotted lines, respec- 
tively. Note the plateau, 
corresponding to the stellar 
"~ 1.6 Mm bump" feat ure 
(ISimpson and Eisenhardtl 
1999). 
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2004). This places them on the upper e nd of the stellar mass functio n over 
the entire redshift range < z < 4 (e.g. iRocca-Volmerange et alll2004T) . 

A ma jor tool in HzRG resea rch is the Hubble K ~ z diagram, intro- 
duced by lLillv and Longain ( 1984T ). Despite significant spectral correction ef- 
fects (K-corrections), radio galaxies form a remarkably narrow sequence out 
to to z ~ 3, especially wh en compared with near-IR selected field galaxies 
(e.g. |Pe Breuck et all l2~002)fFig. [T2]) . The small scatter in the K — z relation 
for HzRGs was found to be weakly correlated with radio luminosity (e.g. 
iBest et alll998| ). Such a correlation can be understood if the radio power is 
dominated by the Eddington limiting lum inosity of the nuclear supermassive 
black hole ( Rawlings and Saunders! fl99lt ) (Section [5.2[) and thus a measure 
of the black hole mass. The well-estab lished relation between black hole and 
bulge mass (e.g. iMagorrian et allll998f ) then implies that radio power would 
be correlated with host galax y mass, as appea rs to be the case for low-redshift 
less- luminous radio galaxies ( Best et a3l2005l) . 

Caution should be exercised in determining the K-band continuum mag- 
nitudes for HzRGs and in interpreting the details of the K-z diagram in 
terms of galaxy masses. There are several reasons for this. First, it is often 
difficult to disentangle the stellar component from the warm dust compo- 
nent of the K-band flux. Secondly, depending on the redshift, there is also 
a contribution to the K-band emission from bright emission lines. Mass and 
age esti mates derived from the unc o rrected colors are therefore oft en overes- 
timated (jEales and Rawlingsl Il996t IRocca-Volmerange et all l2004h . Thirdly, 
the stellar population models on which the masses are based are subject to 
considerable uncertainty. For example, incorporation of recent data on ther- 
mally pulsati ng asymptotic gi ant branch stars reduce the derived masses by 
large factors ( BruzualAll2007l ) . In addition the Initial Mass Function and the 
gas to dust content at early epochs is highly uncertain. 

Because of its ability to measure the rest-frame near-IR fluxes of HzRGs 
radio galaxies, the Spitzer Telescope, h as facilitated c o nsider able progress 
in this field during the last few years. I Seymour et all ( 20071 ) surveyed 69 
radio galaxies having 1 < z < 5.2 at 3.6, 4.5, 5.8, 8.0, 16 and 24 ^m. They 
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Fig. 12 Composite Hubble K-z diagram of radio an d optically selected galaxies, 
slightly modified from Rocca-Volmcrangc ct al (2004) and kindly provided to us by 
Brigitte Rocca. The radio galaxies are denoted by the red squares. The optically- 
selected galaxies are plotted in black (HDFN = crosses and Hawaii survey = dots) . 
The radio galaxies trace the upper envelope of the K-z diagram, with HzRGs being 
amongst the brightest galaxies in the early Universe. Note that continuum fluxes 
derived from K-band or H-band magnitudes can be contaminated by bright emission 
lines in the band. Also plotted are K-magnitude tracks for elliptical galaxies formed 
from various initial reservoir gas masses (baryoni c). These were calculated usin g 
the Pegase galaxy evolution models as described in Rocca- Volmerangc ct al (2004). 

decomposed the rest-frame optical to infrared spectral energy distributions 
into stellar, AGN, and dust components and determined the contribution of 
host galaxy stellar emission at rest-frame H band. By consistently deriving 
the stellar luminosity at the same rest-frame wavelength near the peak of the 
stellar emission at A ros t ~ 1.6 /im, effects of measuring the K — z diagram 
with galaxies over a large range of redshift through fixed observing bands 
were eliminated and contamination of the fluxes by AGN and emission lines 
was minimised. 

The Spitzer results confirm that radio galaxies have stellar masses mostly 
within the 10 11 to 10 12 M range, with almost no dependence on redshift. 
The fraction of emitted light at rest-frame H band from stars was found to be 
> 0.6 for ~ 75% of the HzRGs. The weak correlation of stellar mass with radio 
power determined using ground-based data is only marginally significant in 
the Spitzer results. Rest- frame near-IR studies of less luminous radio sources 
are needed to extend the range in radio luminosity and settle whether there 
is indeed a dependence of stellar mass on radio luminosity. 
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Redshift 

Fig. 13 Rest-frame H-band stellar luminosity vs. re dshift for a sample o f radio 
galaxies observed with the Spitzer Telescope. [From [Seymour et a l (2007)]. The 
masses of stars are derived from the best-fit models to the multiband photometry, 
including Spitzer fluxes. Solid circles indicate HzRGs detected at 24 - 160 /im with 
MPIS. Open circles are luminosities derived from HzRGs that were not detected by 
MIPS. Upper limits are indicated by arrows. Dashed lines correspond to elliptical 
galaxies formed at z/ orm = 10. Crosses mark the stellar luminosity of submillimeter 
galaxies. 



4.2 Young stars - UV absorption lines 



HzRGs undergo vigorous star formation. The most direct evidence for this 
is the detection of rest-frame UV photosphcric stellar absorption lines and 
P-Cygni features driven by stellar winds. Because observation of these lines 
requires several hour exposures on lOm -class telescope s, their detection is 
limited to the S V A1502 line in 4C 41.17 (De v et allll997h. tentativ e measure- 
ment of C III A1428 in TN J2007-1316 (be Breuck and Reulandll2005h and 
VLT spectroscopy of th e Spiderweb Galaxy (jMilev etldl 20061 : iNesvadba et all 
2006; Ha tch et a]||2007f ). The presence of absorption lines are direct evidence 
for prodigious star formation in 4C41.17 (up to 150 M0yr _1 , after correction 
for extinction) (|Dev et alll997l : iReuland et alll2004f ). For radio galaxies whose 
rest-frame UV con tinuum is polarised, scattered AGN continuum emission 
(jVernet et alll200lh . make it difficult to measure the star formation rates 
directly. 

Additional evidence that substantial star formation is occurring in HzRGs 
is provided by the sub- millimeter observations of dust, for whose heating 
star formation rates of thousands Mq per year are needed (but see Section 
I3.5P and the systematically larger Ly a luminositi es and Lyq/Hell ratios fo r 
HzRGs with z > 3, than for those with 2 < z < 3 (jVillar-Martm et all[2007U ) 
(Section [3X2]) ■ 
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Fig. 14 Spectrum of the 
central 2 1 of the HzRG 
4C 41 17 at z = 3 8. 
[From IDev et all (fl997l )]. 
This is compared with the 
UV spectrum of the Bl 
star-forming knot in the 
nearby Wolf-Rayet star- 
burst galaxy NG C 1741 
l |Conti et alfl996t ). The or- 
dinate is labeled with the 
flux density scales for 4C 
41.17 and NGC 1741B1 on 
the left and right axes re- 
spectively. The two spec- 
tra show many similari- 
ties in their absorption line 
properties. 



4.3 The alignment effect - Jet-induced star formation 

One of the most remarkable properties of HzRG is the approximate alignment 
between their radio and optical continuum emissions. This ph enomenon was 



comp l etely unexpected when it was discovered two decades ago ((Chambers et all 
Il987t iMcCarthv et~al Il987t) . The "alignment effect" sets in for redshifts 
z > 0.7. Several mo dels have been p roposed or considered to account for 
this alignments (e.g. lMcCarthvlll993l ). The two most promising models are 
(i) dust scattering of l ight from a hidden quasar at the nucleus of the HzRG 



( Tad hunter et a ll998l) and (ii) trig gering of star formation by the synchrotron 



jet as it propagates outward from the nucleus (|De Younglll989t lReeslll989t 
iBegelman and Ciofnl H989I : iBicknell et all2000h . We have already seen that 
radio jets frequently interact vig orously with the w arm gas in HzRGs. 

Since the extensive review bv lMcCarthvl (11993TI . high resolution imaging 



with the HST has provided important n ew information on the morphologies 
of HzRGs and on the align ment effect. ( Pentericci et all[l999l : IBicknell et all 



l200d : iPentericci eFaill200lh . Images with the HS T have revealed that HzRG 
hosts generally have clumpy optical morphologies (jPentericci et alll9 98. 1999). 

Using HST observations of a sample of 20 HzRGs com bined with similar 
resolution VLA radio maps, IPentericci etai (|l999l l200ll ) have for the first 



time examined the close relationship between the radio and rest-frame UV 
and optical morphologies. They find that the alignment effect extends into 
the rest-frame optical, with the z > 2.5 radio galaxies displaying a much 
more clumpy structure than their low-redshift counterparts, which are well 
represented by de Vaucouleurs prof iles. At a slightly lower r esolution, near-IR 
imaging with the Keck Telescope ( van Breugel et alll998l ) also suggest that 



the morphology changes from aligned clumpy structures at z > 3 to relaxed 
de Vaucouleurs profiles at z < 3. However, the redshift-dependance of lumi- 
nosity and the change in the observing wavelength with redshift complicates 
the interpretation of these results. We note that the IR structures of B2/6C 
galaxies at z ~ 1 appear less aligned with their radio structures than are 
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the brighter 3C galaxies at this redshift, imp lying that the st rength of the 
alignment effect depends on radio luminosity ijEales et allll997t ). 

It is now clear that no single mechanism can completely account for all 
the observed properties of the aligned emission. The occurrence of optical 
polarisation in HzRGs in some but not all HzRGs is evidence that dust scat- 
tering plays a role. However, scattering cannot explain the observed detailed 
optical morphologies and similar bending behaviour that has been observed 
between the optical and radio structures. Neither can dust scattering account 
for the alignment of the continuum emission redwards of the rest frame 4000 
A break. The alignment effect is most likely due to a combination of star 
formation induced by the jet and scattering of hidden quasar light along the 
radio source. 

The most comprehensive stud y of jet-induc e d sta r formation has been 
made for of 4C41.17 at z = 3.8 by Bicknell et all (|2000h . The interaction of a 
high-powered (~ 10 46 erg s _1 ) jet with a dense cloud in the halo of 4C 41.17 
is shown to produce shock-excited emission lines and induce star formation. 
Such shock-initiated star formation could proceed on a timescale of a few 
xlO 6 yr), i.e well within the estimated lifetime of the radio source (few xlO 7 

yr) ■ 

Studies of several low-redshift radio galaxies provide further evidence that 

jet-induced star formation occurs. These include NG C 541 / Minkowski's Ob- 

ject f 3C40) (|van Breugel et allll985l : [Croft et alll2006l). 3C 285 (Ivan Breugel and Devi 
119931 ) and the closest radio galaxy Centaurus A (|Grahamlll998l ). Such effects 
can be expected to be much more prevalent during the era of the Universe 
at z > 2, when every galaxy may well have undergone nuclear activity. 

If the red aligned emission is indeed produced by stars, their ages would be 
> 1 Gyr, implying that this mode of star formation has operated almost since 
the formation of the host galaxies. It would also require that the radio jet was 
active a nd aligned in approxima tely the same position angle for this length 
of time. Ivan Breugel et all ( 19981 ) suggested that the tendency of HzRGs to 
be clumpier at z > 3 is due to the increased importance of jet-induced star 
formation at earlier epochs. We know little about how the first population of 
stars in galaxies were made. Although unconventional, it is tempting to spec- 
ulate that jet-induced star formation played a s ignificant role in producing 
stars during the first 3 Gyr after the Big Bang. iGopal-Krishna et all ( 20041 ) 
have suggested that a large fraction of all proto-galactic material within the 
cosmic web was impacted by the expanding lobes of radio galaxies during 
the quasar era, triggering star formation. 



5 ACTIVE GALACTIC NUCLEUS (AGN) 

"Power exercised in secret, especially under the cloak of national security, is 
doubly dangerous" US Senator, W. Proxmire 

Because the luminous collimated radio emission must have been generated 
in their cores, the nuclei of all HzRGs must have been active at some time 
in their histories. 
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Fig. 15 Evidence fo r jet-induced star formation in 4C41.17 at z — 3.8 [ From 
(Bick nell etai l2pQ0h ]. Montage of three HST images taken th rough the F702W , 
F569W (continuum), and Lya niters. The 8 GHz radio images ofG arilli et all (jl994) 
are superimposed as contours. The absence of significant optical polarisation and 
the details of the optical-radio morphologies cannot be explained by the scattering 
models. The Lya image shows a bright arc-shaped feature near B2 at the apex of 
the edge-brightened UV structure, suggestive of a strong bow shock at the loca- 
tion where th e jet i nteracts with dense ambient gas (see also right hand panels). 
Bick nell et all (|2000l ) suggest that this is manifest through the associated shock- 
excited line emission and star formation in the bifurcated structure. 



5.1 Hidden or dormant quasar - The Unification Church 

The most extreme form of nuclear activity is seen in quasars, for which the 
nuclei are so bright that the galaxies are observed as "quasi-stellar objects". 
There is strong evidence that radio galaxies and quasars are manifestations 
of the same parent population. When we observe objects to be radio galaxies, 
the associated quasars are believed to be either (i) hidden due to our special 
viewing angle and/or (ii) highly variable and observed when they are in a 
dormant state. 

The most widely accepted models for unifying quasars and radio galaxies 
postulate that their observed properties are dependent on the orientation of 
the radio source axes r elative to u s . the o bservers. For a review of unifica- 
tion via orientation see lAntonuccil (|l993h : lUrrv and Padovanl (|l995f ); lUrrvl 



(2004|). One of the key factors in establishing the validity of HzRG unifica- 



tion through orientation has been the signature of a hidden quasar detected 
in som e radio galaxies in polarized light. Polarimetry with the Keck Tele- 
scope ()Vernet et alll20Qll ) (See Section [3.5. 2|) show that HzRGs are usually, 
but not always, highly polarised just redward of Lya and that dust-reflected 
quasar light appears generally to dominate their rest-frame ultraviolet con- 
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tinua. Furthermore, broad emission lines, with equivalent widths similar to 
those seen in quas ars, are detected in the two most polarized HzRGs in the 
IVernet et all ( 200lh sample of nine objects. 

There is a convincing case that an obscured quasar exists in those HzRGs 
for which high fractional polarisations and broad polarised recombination 
lines have been measured. What about HzRGs with low polarisation? For 
these objects, it is possible that active quasars are present that are even 
more highly obscured than the polarised HzRGs. However an equally plausi- 
ble e xplanation is that the quas ars are dormant during the epoch of observa- 
tion ( van Groningen et a 1 1980). We know that quasars are highly variable, 
varying by as much as 2 magnitudes during a day. Quasar var iability has 
been studied on timescales of up to 50-years (|de Vries e t al 2006). The vari- 
ability appears to be intrinsic to the AGNs, and occurs in flares of varying 
time-scales, possibly due to accretion-disk instabilities 

During the lifetime of an extended radio sources (at least a few x 10' 
y), quasar activity can be expected to vary considerably. Knots in the radio 
morphologie s can be int erpreted as tracing bursts of quasar activity in the 
radio nuclei ( Mileyl[r980h . According to such a dormant quasar model, if the 
flare duration is substantially shorter than the quiescent periods the objects 
would be more likely to be observed as radio galaxies than radio-loud quasars. 
Also if the flares recurred on a timescale shorter than the recombination time 
(~ 10 4 - 10 5 y), the excitation of the gas would persist. A recurrently flaring 
quasar could also be a source of excitation for some of the non-radio Lya 
halos (Section T3. 2. 5p . 

In summary, there is strong evidence that quasars exist in the nucleus of 
all HzRGs. Although orientation and varying intrinsic activity both proba- 
bly contribute to the observed properties, their relative importance is still 
unclear. 



5.2 Supermassive black hole - Powerhouse of the AGN 

Extragalactic radio sources, and all active galactic nuclei, are believed to be 
powered by gravitational energy produced by rotat ing supermassive black 
holes (SMBHs) in their nuclei (e.g. lLvnden-Bel]||l969t h Material in the galax- 
ies is accreted onto the SMBHs, converted into kinetic energy and ejected as 
collimated relativistic jets along the rotation axes of the S MBHs to produc e 
the radio sources. For a theoretical review of this field, see 

(BlandfordllQQS)- 



Although there is wide agreement about the general idea, many of the details 
are not well understood. 

Observational evidence for the widespread occurance of SMBHs in galac- 
tic n uclei include (i) spe ctroscopy of the gigamaser in the nucleus of NGC 
4258 ( Mivoshi et allfl995h (ii) spectroscopy and photometry of galactic nuclei 



with the HST leading to the correlati on between the masses of galaxy bulges 

and SMBHs in early type galaxies (e.g. Magorrian et allH998l : iKormendv and Gebhardtl 
l200lfc lHaring and Rbj2004t iNovak et al2006j) , (iii ) the small collim ation scale 
of radio sources, measured with VLBI, (10 17 cm (jjunor et allll999l )) and (iv) 
the similarity of the large-scale and small-scale orientation of some giant 
radio sources, indicating that the 'memory" of collimation axis can persist 
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for a s long as ~ 10 8 y (jFomalont and Milevlll97,4 ISchilizzi et aJll97flt iMilevI 

Because of the interaction of the radio jet with the gas, line widths in 
the central regions of HzRGs cannot be used to derive the masses of the 
presumed SMBHs in their nuclei, as for radio-quiet galaxie s. The empirical 
corre lation of SMBH with the bulges of early type galaxies ( Haring and Rbd 



12004 ) would predict black hole m asses of ~ 10 9 M Q at the centre of a 1O 12 M 
HzRG. iMcLure and Jarvisl (|2002t) have derived similar masses for SMBHs in 



radio-loud quasars from the width of their emission lines. However, such 
estimates should be treated with caution. There is a possibility that the 
widths of the emission lines are widened due to pressure from the radio jets 
Also, little is known about how the SMBHs assemble during the formation 
and evolution of a massive galaxy, such as a HzRG, so extrap olation of results 
from lower redshifts may not be valid. [Silk and Reesl ( 19981 ) have postulated 



that supermassive black holes form within the first sub-galactic structures 
that virialise at high redshift, and are in place before most galactic stars 
have formed. 



5.3 Extinction of the "dinosaurs" 

Luminous quasars and radio galaxies have undergone strong evolution in 
their space densities since z ~ 2 (the peak of the " quasar" era" ) and both 
species are now virtually extinct (Section [L4|). The reason why such extreme 
manifestations of nuclear activity are no longer with us is not fully under- 
stood. One obvious explanation is that the supply of "food" for the SMBHs 
began to diminish as more and mor e of the the gas in massive galaxies was 
converted to stars (e.g. lMencill2004l ). 



6 NATURE OF HzRG HOSTS - MASSIVE FORMING 
GALAXIES 

Viktor Ambartsumian: "It is easy to show applying statistical considerations, 
thta the galaxies NGC 5128 (Centaurus A) and Cygnus A cannot be the re- 
sult of the collisions of two previously independent galaxies " 
Rudolf Minkowski "Results of a detailed investigation of NGC 1275 (Perseus 
A ), reported at this symposium, admit to no other interpretation than that 
this system consists of two colliding galaxies ", I A U Symposium No. 5, Dublin, 
1955 

Having considered the various building blocks of HzRGs in detail, we now 
address the nature of HzRGs and how these objects fit into general schemes 
of galaxy evolution. 

There are several reasons for concluding that HzRGs are massive form- 
ing galaxies. First, we have seen that their large near-IR luminosities imply 
that HzRGs are amongst the most it massive galaxies in the early Universe 
(Section 14. Secondly, HzRG hosts have clumpy UV continuum morpholo- 
gies (Section I4.3[) . as expected from galaxies forming through mergers and 
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in accordance with hierarchical models of massive galax y evolution (|Larsonl 
Il992t iDubinskilHoM iGao et allliaooi ISnringel et aill2005h . Thirdly, the rest- 
frame UV spectra and milli meter SEDs indi cate that HzRGs are undergoing 
copious star formation (e.g. iDev et allfl997l) (Section 14. 2|) . 

Not only are HzRGs massive forming galaxies, but there is also strong 
evidence that they are the progenitors of the most massive galaxies in the 
local Universe, i.e. the giant galaxies that dominate the central regions of rich 
local galaxy clusters. We have seen that HzRGs are generally embedded in 
giant (cD-sized) ionized gas halos (Section [321 In (Section [7]) we shall show 
that they are often surrounded by galaxy overdensities, whose structures have 
sizes of a few Mpc (Section [7]). For all these reasons, HzRGs are the probable 
ancestors of brightest cluster galaxies (BCGs) and cD galaxies. 

What can studying the development of HzRGs tell us about the evolution 
of massive galaxies in general? Recently, two processes have been invoked to 
explain discrepancies between the observed colour and luminosity statistics 
of faint galaxies and the model predictions. In contrast with what might be 
expected from simple hierarchical merging scenarios, the data imply that 
star formation occurs earlier in the history of the most massive galaxies 
than for the less massive ones. This process is called galaxy "downsizing " 
(jCowie et alll996HHeavens et alll 2 0041: 1 Thomas et alll2005l:lPanter et all2007h . 
AGNs are preferentially found in the most massive galaxies. To quench star 
formation in these galaxies, it has been proposed that there is some sort 
of feedback mechanism, in which the AG N influences the star formation 
histo r y of the galaxy. This AGN feedback (|Pi Matteo et all l2005t iBest et all 
120051 : (Springel et alll2005l : ICroton et alll2006l : iHopkins et a]|l2006HBestll20o'7f ) 
provides modelers with additional parameters that, not surprisingly, gives a 
better fit to the data. The potential importance of AGN feedback for galaxy 
formation processes was only realised after it was shown that every galaxy 
may host a supermassive black hole (Section [521), establishing a general link 
between galaxies and AGN. 

As the most massive known galaxies in the early Universe and the seat 
of powerful AGN, HzRGs are important laboratories for studying the pro- 
cesses responsible for downsizing and AGN feedback. In practice, feedback 
scenarios are complicated and involve several physical processes that occur 
simultaneously. We have seen in Sections [2] to [5] that there is considerable ev- 
idence for physical interaction between the various building blocks of HzRGs 
- relativistic plasma, gas, dust, stars and the AGN. As they merge, the satel- 
lite galaxies will interchange gas with the ambient medium in the system. 
The gas will move inwards through cooling flows and accretion and provide 
fuel for the supermassive black hole. The SMBHs generate quasars and rela- 
tivistic jets. Hidden and/or dorma nt quasars heat the dust. The jets to gether 
with superwin ds from starbursts (|Armus et all Il99dl : IZirm et a 1 2005) drive 
gas outwards ( Nesvadba et alll2006l ) and can trigger star formation. Shocks 
will be rampant in the chaotic environments in which these processes are 
competing with each other. 

Two modes have been distinguished in theoretical models for coupling of 
the energy output of the AGN to the surrounding gas, with rather confusing 
names. "Quasar-mode" feedback is defined as the situation when radiative 
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coupling occurs that expels the gas and thereby quenches star form ation 
in th e forming galaxy ( Hopkins et all [200^ 1. In "radio- mode" coupling (Best] 
12007ft . star formation is slowed down due to the mechanical energy of the 
synchrotron radio jets that inhibits cooling of the ambient gas. However, 
such rigid distinction b etween feedback processes are highly simplified (e.g. 
iFu and Stocktonll2007h . 

It is impossible to understand how these complicated AGN feedback pro- 
cesses influences the evolution of massive galaxies from a consideration of 
statistical data alone. 

We shall now illustrate some of these processes and their relevance to 
galaxy formation by considering a specific example of a HzRG in more detail. 



6.1 The Spidcrweb Galaxy - a case study 

The Spiderweb Galaxy (MRC 113 8-262) at a redsh ift of z = 2.2 is one of the 



most intensively studied HzRGs (jMilev et all 120061 ). This object provides a 



useful case study for illuminating several important physical processes that 
may occur generally in the evolution of the most massive galaxies. Because 
the Spiderweb Galaxy is (i) relatively close-by, (ii) one of the brightest known 
HzRGs and (hi) is the HzRG with the deepest HST optical image, it is an 
important laboratory for testing simulations of forming massive galaxies at 
the centers of galaxy clusters. 

This large galaxy has several o f the properties expected for the progen - 



itor of a dominant cluster galaxy (jPentericci et allll997|. |1 998. 2000, 2001). 
The host galaxy is surrounded by a giant Lya halo ([Pentericci et all [20001 : 
Kurk et a]||2004lD and embedded in dense hot ionized gas with an ordered 



magnetic field (jCarilli et all Il998h . The radio galaxy is associated with a 3 
Mpc-sized structure of galaxies, of estimated mass > 2 x 10 14 M , the pre- 
sumed antecedent of a local rich cluster (Section [7]) . 

The beautiful ACS Hubble image of the Spiderweb Galaxy is shown in 
Figure 1161 with the radio source and Lya halo superimposed. The figure 
illustrates the structures of the radio, warm gas and stellar components in a 
relatively nearby HzRG. 

It also provides dramatic evidence that tens of satellite galaxies were 
merging into a massive galaxy, ~ 10 Gyr ago. The morphological complexity 
and dumpiness ag rees qualitatively with predictions of hierarchical galax\ 



and dumpiness agrees qualitatively witn predictions or merarcmcai galaxy 
formation models dLarson 19921: Kauffmann et a| 119931 : iBaugh et all ll998t 
lDubinskilll998l:lGao et alll2004l : ISpringel et alll2005h . and llustratcs this pro- 



cess in unprecedented detail. Lya spectroscopy shows relative velocities of 
several hundred km s — 1 , implying that the satellite galaxies ("flies") will tra- 
verse the 100 kpc extent of the Spiderweb many times in the interval between 
z ~ 2.2 and z ~ 0, consistent with the merger scenario. 

An intriguing aspect of the Spiderweb Gal axy is the presen ce of faint 
diffuse emission between the satellite galaxies ( Hatch et all l2007h . Approxi- 



mately 50% of the ultraviolet light from the Spiderweb Galaxy is in diffuse 
" intergalactic" light, extending over about 60 kpc diameter halo. The lumi- 
nosity in diffuse light implies that the emission is dominated young stars 
with a star formation rate of >80 M@ yr _1 . Under reasonable assumptions, 
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the diffuse emission seen in the Spiderweb Galaxy could evolve into the CD 
envelopes seen in many dominant cluster galaxies seen at low redshifts. 

The total mass of all the flies, in the Spiderweb, derived from their UV 
luminosities (assuming lGyr starbursts), is less t han a tenth of th e mass of 
the whole galaxy obtained from its IR luminosity ( Milev et a]|l2006l h Because 
the UV emission is produced by ongoing star formation and the IR emission 
by old stars, this implies that most of the galaxy mass may already have 
assembled by z ~ 2.2, consistent with downsizing scenarios. 

Merging, downsizing and feedback are all likely to be occurring simulta- 
neously in the Spiderweb Galaxy. Merging is a plausible fueling source for 
the nuclear supermassive black holes that produce the radio sources. Pres- 
sure from these ra dio sources is sufficie nt to expel a large fraction of gas 
from the galaxies (jNesvadba et all 12006). thereby quenching star formation 
( Croton et a jl2006h . Because radio lifetimes are relatively short (few x 10 7 yr), 
all massive ellipticals may have gone through a similar short but crucial radio- 
loud phase during their evolution. 

An unexpected feature of the HST image is that there is a signific ant ex- 
cess of faint satellite galaxies with linear structures (|Milev et afeOOfih . These 
galaxies (linear "flies") have similar morphologies (e.g. chains and tadpoles) 
to the linear galaxies that dominate resolved faint galaxies ('1775 > 24) in 
the H ubble Ultra Deep Field (UDF) (|Elmegreen et all 120051 : IStraughn et all 
120061 ) . Although linear galaxies must be an important constituent of the ear- 
liest galaxy population, their nature is poorly understood. Their presence in 
a merging system is relevant for theories of their formation. In the Spider- 
web Galaxy the motions of the flies with velocities of several hundred km s _1 
through the dense gaseous halo, pe rturbed by superwinds from the nucleus 
( Armus et allll990l ; IZirm et alll2005h and the radio jet, would result in shocks. 
The shocks would then lead to Jeans-unstable clouds, enhanced star forma- 
tion along the direction of motion and to ch ain and tadpole morphologies 
(jTaniguchi and Shiovall200ll : iMilev et a]ll200l . 



7 ENVIRONMENT OF HzRGS - PROTOCLUSTERS 

"We have traced the broad outline of the development in regions of excep- 
tional fruitfulness of the first settled village communities into more populous 
towns and cities" "The city community arose around the altar of the seed 
-time blood sacrifice. " The Outline of History, H. G. Wells, 1920 

Within standard Cold Dark Matter (CDM) scenarios the first stars and 
stellar systems should for m through grav i tation al infall of primordial gas 
in large CDM halos (e.g. IWhite and Reeslll978h . Simulations suggest that 
these halos merge and form web-like networks traced by young g alaxies and 
re-ionized gas (e.g. lBaugh et a]|["l998tlDe Lucia and Blaizotll2007l ). The most 
massive galaxies, and the richest clusters emerge from regions with the largest 
overdensities. 
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Fig. 16 The Spiderweb Galaxy. Deep Hubble image of the core of the MRC 1138- 
262 pr otocluster at z = 2.2 obtained with the Advanced Camera for Surveys. 
[From (|Milev et all [2006)]. Superimposed on the HST image are contours of Lya 
(blue, resolution ~1") obtained with ESO's very Large Telescope (VLT), delineat- 
ing the gaseous nebula and radio 8GHz contours (red, resolution 0.3") obtained 
with NRAO's VLA, delineating the non-thermal radio emission. The gaseous neb- 
ula extends for >200 kpc and is comparable in size with the envelopes of cD galaxies 
in the local Universe. 

7.1 Finding protoclusters using HzRGs 

To investigate the emergence of large-scale structure in the Universe, it is 
important to find and study the most distant rich clusters of galaxies. Con- 
ventional methods for finding distant galaxy clusters usually rely on the 
detection of hot cluster gas using X-ray techniques. These are limited by the 
sensitivity of X-ray telescopes to < 1.3, i.e. much smaller than the redshifts 
of HzRGs. 

Since HzRGs are amongst the most massive galaxies in the early Uni- 
verse, they are likely to inhabit regions that are conducive to the forma- 
tion of rich galaxy clusters. During the mid-nineties observational evidence 
emerged that HzRGs are located in dense cluster-type environments. First, 
the measurement of large radio rotation measures (RM) around some HzRGs, 
indicate d that the host ga l axies are surrounde d by a hot magnetized clus- 
ter gas (|Carilh et all Il997t lAthreva et al 1 19981) (Section IP). Secondly, an 
excess of "companion galaxies" was found near HzRGs dRottgering et all 
ll996tlPascarelle et allll996HLe Fevre et a]|l996tlKnopp and Chambersll997l) . 
These results together with the indications that HzRGs hosts are form- 
ing dominant cluster galaxies (Section [6]) prompted the initiation of direct 
searches for galaxy clusters in the vicinity of HzRGs. 

Conventionally the term "galaxy cluster" refers to a bound structure of 
several hundred galaxies. At z ~ 3, the Universe is only about 2 Gyr old, i.e. 
insufficient time for a galaxy with a velocity of a few hundred km/s to have 
crossed cluster-scale structures. Hence any overdense structure of galaxies 
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observed at high r edshifts must st ill be forming and cannot be virialised and 
bound. Following lOverzie 3 (j2006h . we shall use the term "protocluster" for 
an overdense structure in the early Universe (z > 2), whose properties are 
consistent with it being the ancestor of a local bound galaxy cluster. 

During the last few years there have been several successful direct searches 
for protoclusters around HzRGs, with 8 to lOm-class optical/ infrared tele- 
scopes. One search technique is to use narrow-band imaging to detect emis- 
sion line objects at redshifted wavelengths corresponding to those of the 
target HzRGs. A second technique is to carry out broad-band imaging with 
colours chosen to detect "dropout objects" due to redshifted features in the 
continuum spectral energy distributions. The protocluster candidate galaxies 
detected in these imaging experiments were generally followed up by multi- 
object spectroscopy to confirm their redshifts and their membership of the 
protoclusters. The most important e mission line for establishing the redshifts 
of the protocluster g a laxies is Lva fPentericci et all 19971: iKurk et a 
IVenemans etall l2002i 12004 ICroft et all 120051: IVenemans et all 120051 I2007T) 
Othe r relevant emission lines for such searches for z > 3 are Ha ( Kurk et aJ 
I2004I I) or [OIII]A5007. There are two relevant features in the continuum spec- 
tra of galaxies that are exploited for such protocluster sea rches - the Lyman 
break around the 912A Lyman contin uum discontinuity (|Milev et all 120041 
lOverzier et afeOOoMntema et alll2006l) and the 4000 A break or Balmer brea k 
close to 3648A (|Kurk et alll2004t iKajisawa eTal 120061 : iKodama et alll2007h . 
The various galaxy detection techniques are complementary, because they 
tend to select stellar populations with different ages. For example, Lya ex- 
cess galaxies and Lyman break galaxies have young on-going star-forming 
populations, while the Balmer technique is sensitive to stellar populations 
that are older than a few x 10 8 y. 

The measured overdensities of Lya emitters in the radio-selected proto- 
clusters are facto rs of 3 - 5 larger than the fi eld density of Lya emitters at 
similar redshifts ( Venemans et a 1 l2005ll2007h . A recent study with the HST 
of 4000A break objects in MRC1138-262 at z = 2.2, the protocluster that con- 
tains the Spiderweb Galaxy (Zirm et al. 2007, submitted to Ap. J) shows that 
the overdensity of red galaxies is 6.2, compared with non-protocluster fields. 
The photometric redshifts of galaxies in this field shows a significant " redshift 
spike" for 2 < z < 2.5. Although such observations a t z > 5 are sensitivity 
limited, a significant overde nsity of both Lya em itters ( Venemans et alll2004l ) 
and Lyman break galaxies (|Overzier et a 1 2006) has been established arounc 
TN J0924-2201 at z=5.2, the HzRG with the highest redshift known to date. 

Possible protoclusters have also been detected in non-targeted optical 
surveys, with several overdense regions found at large redshifts in the form 
of "redshift spikes" or "filaments" (e.g. ISteidel et allll998j: iKeel et all 11999k 
Steidel et afeoOOtlMoller and Fvnbdl200lHShimasaku et al 200 



2004 iMatsuda et al 120051 ISteidel et all 120051 : lOuchi et all 1200 



gt lHavashino et all 
M). Although 



caution must be exercised in deriving overdensities from the redshift distri- 
butions alone, because th e effects of peculiar galaxy velocities can influence 
the apparent dumpiness ( Monaco et all 2005). However, the observed struc- 
tures around the HzRGs have two additional ingredient expected of ancestors 
of rich clusters, namely (i) the presence of the HzRG hosts with properties 
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Fig. 17 HzRG-selected pro- 
tocluster around TN J1338- 
1942 at z = 4.1 [Fro m 
IVenemans et~a] ([2003 . 120071 )] . 
Shown are the spatial distri- 
bution of the spectroscopically 
confirmed z = 4.1 Lya emit- 
ters (circles and diamonds), 
additional Lya excess candi- 
dates from narrow-band imag- 
ing (stars) and the radio 
galaxy (square). The circles 
represent emitters with red- 
shift smaller than the me- 
dian redshift and the diamonds 
with those larger than the me- 
dian. The size of the circles is 
scaled according to the relative 
velocity of the object. Larger 
circles and diamonds represent 
a larger relative velocity com- 
pared with the median veloc- 
ity. The structure appears to 
be bound in the northwest of 
the image and unbound in the 
south. 



expected of progenitors of dominant cluster galaxies and (ii) sharp peaks in 
the redshift distributions of Lya emitters that lie close to the predetermined 
redshift of the radio galaxy. Taken together, the measured galaxy overdensi- 
ties combined with the presence of HzRGs is strong evidence that the galaxy 
structures around HzRGs are indeed the ancestors of rich local clusters. 

For some of the HzRG-selected protoclusters there is a dditional evidence 
for a dense environment from exces s counts at millimeter ( Ivison et aj feOOO; 
ISmail et alll2003l IStevens et alll2003l: iDe Breuck et a| 2004t lGreve et all l2007ri 
and X-ray ( Pentericci et al 2002t lOverzier et all 2005T ) wavelengths. Also, the 
X-ray counts indicate an enhanced AGN fraction in protoclusters compared 
to the field. 



7.2 Properties of radio-selected protoclusters 

The largest and most comprehensive HzRG protocluster search project was 
based on a study of 8 radio galaxies, with redshifts ranging from 2.2 to 5.2 

using the narrow-band Lya search technique with the VLT dPentericci et al 

2000; [Kurk et afeOOOHVenemans et all2002[|2004ICroft et alll2005l:IVenemans et all 
200,1 l2007l L" Followup observations have been made of candidates in several 



of the targets using a varie ty of ground-bas e d telescopes as well as the Hub- 
ble a nd Sp i tzer Telescope dKurk et a]||2004ll;jMilev et a]||2004 lOverzier et all 
I2005L 120061 : flntema et alll2006l: iKodama et a]|l2007f) . Typically, ~ 3 Mpc-scale 



regions around the HzRGs were covered in the Lya searches. Lya redshifts 
were determined for 1G8 objects. Almost all these objects are star-forming 
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Fig. 18 HzRG-selected pro- 
tocluster around TN J 1338- 
1942 at z = 4.1 as in Fig- 
urc 1171 [ Fro m rVencman s et all 
(2002, [20071 )]. Velocity distri- 
bution of the confirmed Lya 
emitters. The peak velocity of 
the radio galaxy is indicated 
by an arrow. The solid line rep- 
resents the selection function 
of the narrow-band filter nor- 
malized to the total number of 
confirmed emitters. Note that 
the velocity distribution of the 
detected emitters is substan- 
tially narrower than the filter 
width and centered within 200 
km s _1 of the redshift of the 
radio galaxy. Similar data are 
available for 6 protoclusters. 
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Fig. 19 Morphologies of galaxies in the protocluster TN J 1338- 194 2 at z = 4.1 
with the ACS on the Hubble Space Telescope [From [Qverzier et all p006t )]. The 
filter is 1775 and each image measures 3" x 3", corresponding to ~ 20 kpc x 20 kpc 
at z ~ 4. Left are images for a selection of the overdense (7475-dropout objects that 
are candidate Lyman break galaxies (LBGs) in the protocluster. Right are images 
of 12 spectroscopically confirmed Lya excesses (LAE) protocluster members. There 
are a wide range of morphologies present, with some galaxies showing clear evidence 
for small-scale interactions. 



galaxies at similar redshift to the HzRG, with typical star formation rates of a 
few Mq per year, derived from the UV continuum and the Lya luminosities. 

Six of the 8 fields were found to be overdense in Lya emitters by a factor 3 
- 5 as compared to the field density of Lya emitters at similar redshifts. The 
6 included all of the most radio- luminous objects (i.e . all HRzGs with radi o 
luminosities L2.7GHZ > 6 x 10 33 erg s _1 Hz -1 sr _1 ) dVenemans et alll2007lL 
Because these targets were chosen arbitrarily, (redshifts that place Lya in an 
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available narrow-band filter), it is possible that all HzRGs with such large 
radio luminosities are embedded in such overdensities. 

Measured velocity dispersions are in the range ~ 300 - 1000 km s _1 , 
centered within a few hundred km/s of the mean velocity of the radio galaxies. 
Taking account of the relatively narrow peaks in the velocity distributions 
of the Lya emitters compared with the widths of the imaging filters, the 
measured galaxy overdensities range from ~ 5 - 15. Estimates for the sizes of 
the protocluster structures arc limited by the typical size of the imaging fields 
used in the searches to date (~ 8' equivalent to ~ 3 Mpc). Nev ertheless, the 
protocluster sizes are estim ated to be in the range ~ 2 to 5 Mpc ( Interna et all 
120061: IVenemans et al|[2007h . 

Because cluster-size overdense structures at high-redshifts cannot be old 
enough to have become bound, the usual method of calculating cluster masses 
using the virial theorem cannot be applied. However, estimates for the masses 
of the protoclusters can be obtained from the volume occupied by the over- 
density, the mean density of the Universe at the redshift of the protocluster, 
the measured galaxy overd ensity, and the bias parameter (e.g. ISteidel et all 
119981: IVenemans et"a3|2007l\ . The masses obtained a few times 10 14 Mq - 
10 15 M . (iKurk et al 12004 IVenemans et al 120051) ) are comparable to the 
masses of local rich clusters. 

It is interesting to inquire what is the relation of the protocluster struc- 
tures to the general large-scale structure of the Universe. Is the topology of 
the protocluster filamentary, or do HzRG-selected protoclusters illuminate 
the densest most tangled regions of the cosmic web. For TN J1338-1942 at 
redshift z = 4.1, there are indications that the Mpc-sized protocluster of 
Lya excess galaxies may be part of a larger-scale structure. A 25' x 25' sur- 
vey for Lyman break galaxy candidates (B-band dr opouts) showed se veral 
significant density en hancements amidst large voids ([Interna et alll2006i ) and 
lOverzier et a 3 (j2006j)showed that this large scale structure ties in closely with 
significant sub-clustering across the 3' x 3' ACS field near the radio galaxy. 

Future wide-field narrow-band imaging and spectroscopy around HzRG- 
selected clusters should allow the topology of the Universe in the region 
of protoclusters to be mapped in detail, providing a glimpse of large-scale 
structure emerging in the early Universe for comparison with computer sim- 
ulations. 



7.3 Are radio-selected protoclusters typical? 

We have seen that protocluster-like structures have been found around al- 
most all of the most luminous radio galaxies with z > 2 that have been 
targeted. However only ~ 50% of powerful radio sources at z ~ 0.5 are lo- 
cated in rich clusters and radio sources appear to avoid clusters at low redshift 
(|Hill and Lillvlll99lh . What is known about the environment of radio galaxies 
having redshifts 1 < z < 2? This range is difficult to study spectroscopically 
because of the redshift desert (Section II. 3|) . However, there have been sev- 
eral reports of detections of clust ers and/ or excesses of red galaxies around 
radio AGN in this redshift range (|Hall et alll 19981 iKodama and Bowerll2003t 
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iNakata et a] 120021 : iBest et alll2003h . We note that there is a strong correla- 
tion between the radio luminosities and redshift of objects in such studies 
due to Malmquist bias. Taking all the data together, there appears to be a 
substantial increase in the density of the environment around radio galaxies 
as a function of redshift and/or luminosity. 

The statistics of the luminosity function of radio galaxies are consistent 
with every brightest cluster gal axy (BCG) having gone thr ough a luminous 
radio phase during its evolution (|Venemans et al|[2002l 120071 ) . This statement 
is based on the facts (i) that the space density of luminous steep-spectrum 
radio sources decreases by ~ 100 between 2.5 > z > (Section fl~4|) and (ii) 
that the radio synchrotron lifetimes (few xl0 7 yr) are ~ 100 smaller than 
the cosmological time interval corresponding to the observed redshift range 
(few xl0 9 yr). Distant radio galaxies may therefore be typical progenitors of 
galaxies that dominate the cores of local clusters. Likewise, the ancestor of 
every rich cluster in the local Universe may have gone through a phase in 
which it hosted a HzRG. This would imply that HzRG-selected protoclusters 
are typical ancestors of local galaxy clusters. 



7.4 Protocluster evolution 

Radio-selected protoclusters are powerful laboratories for tracing the emer- 
gence of large scale structure and for studying the evolution of galaxies in 
dense cluster environments. An interesting parameter that can be measured 
for line-emitting galaxies is the velocity dispersion of the protoclusters. Al- 
though more statistics ar e needed, the velocity d ispersion appears to decrease 
with increasing redshift ( Venemans et alll2007t) . consist ent with the p redic- 



tions from simulations of forming massive clusters (e.g. lEke et all ll996). 

A population study of the protoclus ter around t he Sp iderweb Galaxy PKS 
1138-262 at z = 2.2 was carried out bv lKurk et all (|2004l l), using deep optical 
and infrared observations with the VLT. Besides Lya emitters, the study 
included objects having apparent Ha excesses and extremely red objects, 
with colours characteristic of old galaxies with 4000 A breaks at the redshift 
of the protocluster. An intriguing result of this study is that candidate Ha 
emitters and 4000 A break objects appear more concentrated towards the 
centre of the protocluster than the Lya emitters. This indicates that the 
galaxies that are dominated by old stellar populations are more settled into 
the gravitational potential well of the protocluster, consistent with them 
being older. 

An important diagnostic of galaxy and cluster evolution is the cluster 
colour-magnitude diagram (CMD). The presence of a narrow 'red sequence" 
in the cluster CMD is well established out to a redshift of 1.4 ( Sta nford et all 



2005). The red colours imply that the galaxies are dominated by older stars, 
whose SED peaks redwards of the 4000A break. The presence of this red 
sequence can be used to set a lower limit to the redshift at which the stellar 
populations formed. It is extremely difficult to measure the CMD at z > 1.4, 
even with the largest ground-based telescopes. 

However, indications of an emergent red sequence has been found in the 
CMD of the protocluster MRC1138-262 at z = 2.2, both with a wide-field 
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near-IR imager on the Subaru Telescope ( Kodama et all (2007); Fig. l20l and 
with the much narrow field but more sensitive NICMOS on the HST (Zirm 
et al. 2007). The galaxy colours indicate that while some relatively quiescent 
galaxies exist in the prot ocluster, most of the g alaxies are still undergoing star 
formation Furthermore. iKodama et~al ( 2007 ') studied 3 other radio-selected 
protoclusters and found that the fraction of red galaxies in the 4 protoclusters 
increases between z ~ 3 and z ~ 2. To summarise, there is strong evidence 
that most protocluster galaxies undergo substantial star formation between 
z ~ 3 and z ~ 2 and that the bright end of the red sequence is still being 
formed during this epoch. 

Information about the early evolution of stellar populations in clusters has 
been obtained at even higher redshifts. In the protocluster surrounding the 
HzRG TN J1338-1942 at z = 4.1, the distribution of Lya emitters (LAEs) is 
highly fila mentary and appears to avoid the locations of Lyman break galaxy 
candidates ( Overzier et alll20"06T) . A similar spatial segregation between LAEs 
and LB Gs was observed in a structure around QSO SDSS J0211-0009 at 
z=4.87 (jKashikawa et a This indicates that an age- or mass-density 

relation was emerging little more than 1 Gyr after the Big Bang, when the 
Universe was only 10% of its present age. 

We note that derived masses of z > 2 LBGs and LAEs are <; 10 smaller 
than the masses of early-type galaxies in local clusters, indicating that a 
large fraction of the stellar mass still has to accumulate through merging 
( Overzier et alll2b06t ). Detailed observations of protocluster regions on much 
larger scales (~50 co-moving Mpc) are needed to test if the number density of 
LBGs is indeed consistent with forming the cluster red sequence population 
through merging. Simulations show that clusters with masses of > 10 14 M 
can be traced back to regions at z —4 - 5 of 20 - 40 Mpc in size, and that 
these regions are associated with overdensities of typical dark matter halos 
hosti ng LAEs and LBG s of 8 g ~ 3 and mass overdensities 6 m in the range 0.2 
- 0.6 ( Suwa et""aJ 12006). Also, recent numerical simulations of CDM growth 
predict that quasars at z ~ 6 may lie in the center of very ma ssive dark 
matter halos of - 4 x 10 12 M Q ([Springel et alll2005tlLi et alll2007D . They are 
surrounded by many fainter galaxies, that will evolve into massive clusters 
of - 4 x 10 15 M Q at z = 0. 

The addition of the new IR - oprical camera WFC3 to the HST and the 
advent of several new wide-field imagers and multi-object spectrographs on 
ground-based telescopes will allow the detailed evolution of protoclusters to 
be studied in great detail during the next few years. 



8 CONCLUSIONS AND FUTURE 

"Curioser and curiouser! cried Alice (she was so much surprised that for the 
moment she quite forgot how to speak good English). Now I'm opening out 
like the largest telescope that ever was!" Alice's Adventures in Wonderland, 
Lewis Caroll, 1865 

In this review we have shown that HzRGs are fascinating objects in their 
own right and that they provide important diagnostics for studying the early 
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Fig. 20 Color-magnitude diagram o f J - Kg plotted again st Ks for the MRC 1138- 
262 protocluster at z — 2.2. [From Koda ma et all ([2007)]. Filled circles indicate 
protocluster member candidates selected using colour selection criteria. Large stars 
mark the targeted radio galaxies. Dotted error bars show Is photometric errors. 
Solid, dashed and dotted lines show the expected location of the CMD at the 
relevant redshift for passive evolution. The iso-stellar mass lines of lO n M0 and 
10 10 Mq are shown by thick dotted lines. The dot-dashed lines indicate 5 a (Kg) 
and 2cr (J) detection limits. 



Universe. Knowledge about HzRGs and their evolution is fundamental to 
understanding the formation and evolution of galaxies and the large scale 
structure of the Universe. 

There are still many aspects of HzRGs that are not understood. Here are 
some open questions. 

— What is the particle acceleration mechanism in the relativistic plasma 
and why do HzRGs have much steeper radio spectra than nearby radio 
sources ( Section HOI) . 

— Why caused the most luminous radio galaxies and quasars to become 
virtually extinct between z ~ 2 and the present (Sections 11.41 and !5.3p ? 

— What is the nature, extent and kinematics of the hot gas, one of the most 
massive and least studied constituents of HzRGs (Section 13. ip ? 

— What is the origin of the Lya halo that appears to be falling into the 
HzRG and how is this warm gas (filling factor ~ 10~ 5 ) distributed with 
respect to the hotter and colder gas (filling factors ~ 1) (Sections 13.2.41 
and El? 
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— What are the detailed processes by which the radio jets interact with the 
gas and trigger starbursts and how important is jet-induced star formation 
for producing stars in the early Universe (Sections 13.2.31 and !4.3|) ? 

— What are the temperature and densities and composition of the molec- 
ular gas and dust? how is this distributed spatially and what are the 
implications for the star formation histories (Sections 13.41 and 13. 5ft ? 

— What is the relative importance of AGN and starburst heating in the 
mid-IR emission of HzRG? 

— What effect does galaxy merging have on star formation and what phys- 
ical effects are responsible for downsizing (Section [6])? 

— What effect does feedback between the AGN and the galaxy have on the 
evolution of HzRGs and the general evolution of massive galaxies (Section 

— What is the detailed mechanism by which the SMBHs produce quasars 
and luminous collimated jets? 

— How do the AGN/quasars vary on long time-scales and what are the 
causes of the variability (Sections 15. II and 15. 2j) ? 

— How is the SMBH built up and what role does merging play in this 
evolution fSection !5.2p ? 

— What is the size distribution of radio-selected protoclusters and what is 
is the topology of the cosmic web in the neighbourhood of the HzRGS 
(Section El)? 

— How do the various populations and constituents of distant protoclusters 
evolve and eventually become virialised clusters (Section 17. 4p ? 

— What is the radio luminosity function of galaxies in radio-selected proto- 
clusters? Is the radio emission from the protocluster galaxies an impor- 
tant contributor to forming the radio halos, that are seen at the centres 
of many nearby rich clusters? What is the strength and configuration of 
the protocluster magnetic fields (Sections [7] and [2])? 

There are good prospects for making progress during the next decade. 
Since its inception, the study of radio galaxies, has been observationally 
driven. Several forefront astronomical facilities now being constructed or 
planned that will give new insights into the nature of HzRGs and their envi- 
ronments. 

First, with a combination of sensi tivity and spatial reso lution, the new 
low-frequency radio arrays, LOFAR ( Rottgering et all 120061 ) and the LWA 



(|Kassim et all [20061. will open up the frequency window below ~ 50 MHz 



for HzRG studies. LOFAR will survey the sky to unprecedented depth at 
low-frequencies and will therefore be sensitive to the relatively rare radio 
sources that have extremely steep spectra. Because of the a vs z relation 
(Section EH), LOFAR is likely to detect HzRGs at z ~ 8, if they exist. Stud- 
ies of detailed low-frequency spectra and their spatial variations will provide 
new information about the mechanism responsible for the a vs z relation. 
Presently combination of the new radio surveys with planned new deep op- 
tical and infrared wide-field surveys, such as PAN-STARRS (iHodapp et all 
l2004h and those with the VST and VISTA (lArnaboldi et~al l2007h will be 



used to identify HzRGs and provide photometric data. 
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Another task for sensitive radio arrays, such as LOFAR, the EVLA (|Napierl 
120061 ) and eventually the Square Kilometre Array (SKA), will be to survey 
the radio emission of galaxies in protoclusters. The new arrays will study 
radio emission produced by relativistic jets and be able to detect and inves- 
tigate radio emission from the brightest star forming protocluster galaxies. 

Secondly, ALMA (jBrown et a]ll2004l ) and the EVLA, with their unprece- 
dented sensitivities and resolutions at millimetre and sub-millimetre wave- 
lengths, will revolutionise the study of molecular gas and dust. Several dif- 
ferent CO transitions can be observed, allowing entire "CO ladders" to be 
constructed and the density and temperature structure of the molecular gas 
to be unraveled. Fainter molecular lines can be used to trace even denser gas 
than that studied until now. Important information about the dust compo- 
sition and the gas to dust ratios is likely to be obtained. 

ALMA's sensitivity at millimetre wavelengths should also facilitate ob- 
servations of the atomic CI Carbon lines in HzRGs. This would provide an 
important constraint on the global metalicity of the gas. The fine-structure 
line of C + at A res t=157.74/im line is one of the main cooling lines in nearby 
galaxies, and has now also been detected in several of the most distant quasars 
known. 

Thirdly, we can expect considerable progress in disentangling the detailed 
evolutionary history both of HzRGS and of radio-selected protoclusters. This 
evolutionary detective work will be pursued by combining spectroscopic data 
from the next generation of spectrograp hs on 8m-class telesc opes with imag- 
ing results from the new camera, WFC3 ([Leckrone et al|[T998l ). on the Hubble 
Space Telescope. For example, the detection of supernovae in z ~ 2 protoclus- 
ters will become possible. On a longer timescale, tracing the detailed history 
of the formation and evolution of HzRGS and the surrounding protoclusters 
will be helped enormously with the advent of 30m -class ground-based tele- 
scopes in the optical and near -infrared and the James Webb Space Telescope 
(JWST) (ICardner et all200fih in the near and mid-infrared. 

F ourthly, the next gene ration X-ray telesc ope, such as XEUS ( Bavdaz et all 
120061 ) or Constellation-X (|Tananbauml [20 06:) , will make observations of X- 
rays from HzRGs an important tool for studying galaxy formation. It will 
have sufficient sensitivity to perform spectroscopic studies of hot gas in 
HzRGs. The 0.3 - 10 keV X-ray band contains the inner (K-shell) lines for 
all of the abundant metals from carbon to zinc as well as many L-shell lines. 
These atomic transitions provide important new plasma diagnostics of the 
HzRG hot gas. 

Fifthly, and perhaps most exciting, the potential discovery of HzRGS with 
z > 6 could open up a unique new window for studying the very early Uni- 
verse during the " Epoch of Reionisation" . Recent observational constraints 
on suggest that cosmic reio nization may have taken place between z ~ 11 
and z ~ 6 (|Fan et aJl20"53 ). The existence of HzRGs within the near edge 
of cosmic reionization could be used as sensitive probes of intermediate- to 
small-scale structures in the neutral IGM, through redshifted HI absorptio n 
observations (Section 13.3. 1|) (|Car illi et all l2002t iFurlanetto and Loebl [2002) , 
complementary to the very large scale that can be studied in HI emission. 
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Finally, we point out that in trying to understand HzRGs or in using them 
to probe the properties of the early Universe, we should exercise humility, 
remembering that 96% of the energy density in the universe is in a form of 
dark matter and dark energy, that are not directly observable. 

HzRGS pinpoint the progenitors of cosmic megacities in the local Uni- 
verse. Studying them is akin to investigating how human civilisations were 
urbanised. During the last few decades we have obtained some glimmerings 
of how and why this cosmic urbanisation occurred. During the next few 
decades we shall learn a great deal more about the origin of the megacities, 
the nature of the first city dwellers, their organisational structure and their 
successes and failures. 
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9 APPENDIX: KNOWN HZRGs 

This appendix is a compendium of HzRGs (z > 2, Li5oo(rest) > 10 27 ' 5 erg 
s _1 Hz -1 ) known at the time of writing. Table 3 lists the HzRGs, their 
positions, redshifts, radio luminosities extrapolated to 500MHz rest frequency 
and relevant references. For convenience the HzRGs are listed in order of right 
ascencion in Table 4. 

The luminosities were calculated using 

L v = AnDlS^l + z)- a ~\ 

where Y>l is the luminosity distance Q, S v is the observed flux density at the 
rest frequency, and a is the spectral index (S v oc v a ). 

Why did we choose a rest-frame frequency of 500 MHz to express the 
radio luminosities? The radio sources have steep spectra and are generally 
taken from low- frequency surveys, so it is appropri ate to express their lu- 
mino sities at as low a rest frequency as possible (e.g. iBlundell and Rawlingsl 
l2000h . For consistency, surveys covering the northern and southern sky with 
comparable spatial resolution should be used. The most a ppropriate surveys 
are the 74 MHz VLA Lo w-frequency Sky Su rvey fVLSS; ICohen et alll2007t) 
and the 1.4 GHz NVSS (jCondon et all 1 19981) . These contain ~ 75% of the 
known HzRGs. For the entire redshift range of the sample, the observed fre- 
quency associated to i/ rcst =500MHz is located close to the VLSS frequency, 
but still located between VLSS and NVSS, thereby avoiding extrapolations 
to un-observed parts of the radio spectrum. 

For HzRGs not included in these catalogues, th e following combinati on of 
surveys were u sed in order of prefe rence: WENSS (iRengelink et all 1997) and 



NVSS, Texas (iDouglas et allll996l) an d NVSS, S UMSS (iMauchet all 12 003) 
and NVSS, or MRC ([Large et allll98l[ ) and PMN (jGriffith and Wrightjll993¥ . 

1 assuming ffo=70kms _1 Mpc -1 , iXi=0.7 and J?m=0.3 
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A small number of HzRGs from the CENSORS survey (jBrookes et a]|l2006h 
do not yet have measured radio spectra. For these we assumed a constant 
spectral index of a = -1.1, leading to uncertainties in the derived luminosities 
that could be as much as 0.5 dex. The relevant luminosities are indicated by 
the daggers. 

The table also lists references for the redshifts and relevant followup ob- 
servations. Although references to the followup work are not complete, they 
provide a useful pointer for readers who require further information. The key 
to the references to followup observations are: CO = CO, G = General, H 
= HST imaging, HI = HI absorption, M = Millimetre/ submillimetre, N 
= Nebular continuum, OP = Optical polarisation, OS = Optical/infrared 
spectroscopy, P = Protocluster, R = Radio, S Spitzer IR and X = X-rays. 



Table 3 

Known radio galaxies with z > 2 and radio luminosities Lsoofre 



RA(.T2000) DEC(J2000) 



LogtLsooT" 



10 27 - 5 W RV 



Name 



MRC 2036-254 

MRC 0015-229 

MRC 0156-252 

MG 2348+0507 

MP J1758-6738 

TXS 0448+091 

NVSS J002402-325253 

MRC 2048-272 

B3 1204+401 

PKS 0011-023 

5C 7.223 

MG 1401+0921 

7C 1740+6640 

7C 1802+6456 

TN J1102-1651 

4C +39.21 

6C 1134+369 

TXS 2034+027 

TXS 0214+183 

WN J1242+3915 

4C +40.49 

NVSS J103615-321659 
4C +43.31 

NVSS ,1144932-385657 
TXS 0355-037 
TNR 2254+1857 
MRC 1138-262 

6C**0746+5445 
6C* 0024+356 
MRC 0030-219 
B2 1056+39 
6C* 0135+313 
5C 7.10 
5C 7.269 
5C 7.271 



20:39:24.5 
00:17:58.2 
01:58:33.5 
23:48:32.0 
17:58:52.9 
04:51:14.6 
00:24:02.3 
20:51:03.4 
12:07:06.2 
00:14:25.5 
08:24:06.4 
14:01:18.3 
17:40:42.2 
18:02:22.1 
11:02:47.1 
07:58:08.8 
11:37:07.7 
20:36:34.8 
02:17:25.8 
12:42:53.1 
23:07:53.5 
10:36:15.3 
13:52:28.5 
14:49:32.8 
03:57:48.0 
22:54:53.7 
11:40:48.3 

07:50:24.6 
00:26:52.0 
00:33:23.9 
10:59:11.5 
01:38:06.6 
08:11:26.5 
08:28:39.7 
08:28:59.5 



-25:14:31 
-22:38:04 
-24:59:30 
+05:07:33 
-67:38:34 
+09:14:29 
-32:52:54 
-27:03:05 
+39:54:39 
-02:05:56 
+26:28:12 
+09:21:24 
+66:38:44 
+64:56:45 
-16:51:34 
+39:29:28 
+36:39:54 
+02:56:55 
+ 18:37:03 
+39:15:49 
+40:41:49 
-32:16:57 
+42:59:23 
-38:56:58 
-03:34:08 
+ 18:57:04 
-26:29:10 

+54:38:07 
+35:56:24 
-21:42:01 
+39:25:01 
+31:32:42 
+26:18:19 
+25:27:30 
+24:54:00 



2.000 

2.01 

2.016 

2.019 

2.026 

2.037 

2.043 

2.06 

2.066 

2.080 

2.092 

2.093 

2.10 

2.110 

2.111 

2.119 

2.125 

2.129 

2.131 

2.131 

2.140 

2.144 

2.149 

2.152 

2.153 

2.154 

2.156 

2.156 
2.161 
2.168 
2.171 
2.199 
2.185 
2.218 
2.224 



28.52 
28.32 
28.46 
28.58 
28.99 
28.11 
27.93 
28.72 
28.33 
28.15 
27.71 
28.26 
27.77 
28.18 
28.07 
28.51 
28.31 
28.45 
28.46 
28.03 
28.46 
27.82 
28.43 
28.07 
28.34 
27.80 
29.07 

27.64 
28.04 
28.14 
28.20 
27.97 
28.21 
27.78 
27.97 



References 

■R^rorsriopx 2301 - 

RI43I: m n| \10 



hEIR-EI? ;HE3? ] S[H1 X [230l 



R[246TlSf280llXf230l 



R 246] 
R I2461 

OS 320]R 43] 



os cm 



G[2llH[24^[^ :[2T lMp9S[263t OS f22l P [2451 [T80l : [2481 : \wj [l8ll : 

63; 176; 325]R 43; 9; 8; 10]S 280]X 44; 46] 



Ml 



S2 Spinrad et al MG list (private communication) 



Name 


RA(J2000) 


DEC(J 


2000) 


z 


Ref. 


6C* 0142+427 


01 


45 


2!) 





+42 


57 


42 


2 


225 


157 


TXS 0200+015 


02 


02 


42 


9 


+01 


49 


11 


2 


229 


271 


NVSS J101008-383629 


10 


10 


08 





-38 


36 


29 


2 


236 


38] 


TXS 1113-178 


11 


16 


14 


7 


-18 


06 


23 


2 


239 


271 


6C* 0115+394 


01 


17 


55 


3 


+39 


44 


33 


2 


241 


157 


6C 0629+53 


06 


33 


52 





+53 


16 


31 


2 


246 


102 


6C 0901+54 


09 


05 


24 


2 


+54 


05 


42 


2 


249 


102 


TN J0452-1737 


04 


52 


26 


7 


-17 


37 


54 


2 


26 


72] 


4C 40.36 


18 


10 


55 


7 


+40 


45 


23 


2 


265 


48], 






















327 


6C* 0152+463 


01 


55 


45 


6 


+46 


37 


11 


2 


279 


157 


4C 40.22 


08 


59 


59 


6 


+40 


24 


36 


2 


28 


234 


MRC 1324-262 


13 


26 


54 


7 


-26 


31 


43 


2 


28 


206 


MG 1747+1821 


17 


47 


07 





+ 18 


21 


11 


2 


281 




6C* 0106+397 


01 


09 


25 


4 


+40 


00 


01 


2 


284 


157 


MP J0340-6507 


03 


40 


44 


6 


-65 


07 


12 


2 


289 


32] 


6C 1106+380 


11 


09 


28 


9 


+37 


44 


31 


2 


290 


102 


MG 0936+0503 


09 


36 


10 


1 


+05 


03 


49 


2 


306? S2 


MG 1251+1104 


12 


51 


00 


1 


+ 11 


04 


20 


2 


322 


298 


MRC 0349-211 


03 


51 


11 


8 


-20 


58 


00 


2 


329 


201 


TXS 0211-122 


02 


14 


17 


4 


— 11 


58 


47 


2 


340 


318 






















327 


4C 48.48 


19 


33 


05 


6 


+48 


11 


46 


2 


343 


51] 


TXS 1707+105 


17 


10 


06 


f) 


+ 10 


31 


09 


2 


349 


271 


BRL 0128-264 


01 


30 


27 


9 


-26 


09 


57 


2 


348 


18] 


6C* 0118+486 


01 


21 


16 


3 


+48 


57 


40 


2 


350 


157 


NVSS J094949-213432 


09 


49 


49 





-21 


34 


33 


2 


354 


35] 


PKS 1425-148 


14 


28 


41 


7 


-15 


02 


28 


2 


355 


72] 


NVSS J232651-370909 


23 


26 


51 


5 


-37 


09 


11 


2 


357 


79] 


4C -00.54 


14 


13 


15 


1 


-00 


23 


00 


2 


360 


271 


MG 0001+0846 


00 


01 


15 


5 


+08 


46 


39 


2 


36? 


327 
32 


B2 1159+39A 


12 


01 


50 





+39 


19 


11 


2 


370 


114 


6C**0854+3500 


08 


57 


15 


9 


+34 


48 


24 


2 


382 


66] 


LBDS 53W002 


17 


14 


14 


8 


+50 


15 


30 


2 


390 


351 


6C 0930+389 


09 


33 


06 


9 


+38 


41 


48 


2 


395 


102 



' Radio luminosity calculated assuming a spectral index a = 
S2 Spinrad et al MG list (private communication) 
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Table 3, 



Name 


RA(J2000) 


DEC(,T2000) 


z I 


4C 34.34 


11 


16 


30 


4 


+34 


42 


24 


2 


400 


7C 1736+650 


17 


36 


37 


5 


+65 


02 


28 


2 


400 


TXS 0748+134 


07 


51 


01 


1 


+ 13 


19 


27 


2 


419 


NVSS J095226-200105 


09 


52 


26 


5 


-20 


01 


05 


2 


421 


TN J1033-1339 


10 


33 


10 


7 


-13 


30 


52 


2 


425 


NVSS J094925-203724 


09 


49 


25 


9 


-20 


37 


24 


2 


427 


4C 40.02 


00 


30 


49 





+41 


10 


49 


2 


428 


B3 0731+438 


07 


35 


21 


9 


^43 


44 


21 


2 


429 


MRC 1106-258 


11 


08 


30 


3 


— 26 


05 


06 


2 


43 


5C 7.15 


08 


12 


19 


2 


+26 


30 


18 


2 


433 


MRC 0406-244 


04 


08 


51 


4 


—24 


18 


17 


2 


44 


6C**0834+4129 


08 


37 


49 


2 


+41 


19 


54 


2 


442 


8C 1536+620 


15 


37 


09 


5 


^61 


55 


36 


2 


450 ] 


MG 2308+0336 


23 


08 


25 


o 


+03 


37 


04 


2 


457 


4C 12.32 


08 


55 


21 


4 


+ 12 


17 


27 


2 


468 


3C 257 


11 


23 


09 


4 


+05 


30 


18 


2 


474 


4C 23.56 


21 


07 


14 


g 


+23 


31 


45 


2 


479 


TXS 2332+154 


23 


54 


58 


(} 


+ 15 


45 


50 


2 


481 


MRC 2104-242 


21 


06 


58 


2 


—24 


05 


11 


2 


491 


WN J0303+3733 


03 


03 


26 




^37 


33 


42 


2 


505 


TXS 1647+100 


16 


50 


05 


f) 


+09 


55 


11 


2 


509 


NVSS J213510-333703 


21 


35 


10 


5 


-33 


37 


04 


2 


518 


TXS 1558-003 


16 


01 


17 


4 


-00 


28 


46 


2 


520 


8C 1039+681 


10 


42 


37 


4 


+67 


50 


24 


2 


530 


NVSS J231144-362215 


23 


11 


45 


2 


-36 


22 


15 


2 


531 


6C* 0112+372 


01 


14 


50 


2 


+37 


32 


33 


2 


535 


TXS 1436+157 


14 


39 


04 


9 


+ 15 


31 


19 


2 


538 


WN J1115+5016 


11 


15 


06 


9 


+50 


16 


24 


2 


54 


MRC 2139-292 


21 


42 


16 


7 


-28 


58 


38 


2 


55 


MP J1755-6916 


17 


55 


29 


8 


-69 


16 


54 


2 


551 


TXS 2319+223 


23 


21 


42 


4 


+22 


37 


55 


2 


554 


6C**1045+4459 


10 


48 


32 


4 


+44 


44 


28 


2 


571 


TXS 0828+193 


08 


30 


53 


4 


+ 19 


13 


16 


2 


572 



Rcf. 



LR Lacy and Rawlings (private communication) 
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Table 3, continued 



Name 


RA(J2000) 


DEC(J2000) 


z 


Ref. 


L°g(^5aoMi 


z References 


PKS 0529-549 


05:30:25.4 


-54:54:21 


2 


575 


271 


29.16 


OSI32CIR.I34IISI280II 


TXS 2353-003 


23:55:37.0 


-00:02:58 


2 


592 


34] 
7j 


28.76 




WNH 1736+6504 


17:36:37.4 


+65:02:35 


2 


6 


25al 


27.74 




MRC 2052-253 


20:55:01.5 


-25:09:53 


2 


60 


2091 


28.77 


Rl24l 


WN J0040+3857 


00:40:56.2 


+38:57:30 


2 


606 


3T 


27.80 




4C 26.38 


12:32:23.6 


+26:04:07 


2 


609 


51] 


28.41 


B3 0811+391 


08:14:30.6 


+38:58:35 


2 


621 


233] 


27.84 




PKS 0742+10 


07:45:33.0 


+ 10:11:13 


2 


624 


m 


28.74 




MRC 2025-218 


20:27:59.5 


-21:40:57 


2 


63 


201 


28.74 


H[243;[243 OS[34lR[4l|i;[|;[ll Sf280l 


6C* 0122+426 


01:25:52.8 


+42:51:51 


2 


635 


1571 


28.04 




MRC 0140-257 


01:42:41.2 


-25:30:34 


2 


64 


206] 


28.44 


Hf247lOS[34lR[43llgtll[ulS 


TN J1941-1952 


19:41:00.1 


-19:52:14 


2 


667 


3j 


28.85 




PKS 1357+007 


14:00:21.4 


+00:30:19 


2 


671 


27J 


28.72 


OSf32GlR[24l 


NVSS J233034-330009 


23:30:34.5 


-33:00:12 


2 


675 


Hp 


28.23 




7C 1758+6719 


17:58:47.2 


+67:19:46 


2 


70 


1841 


27.99 




TXS 2202+128 


22:05:14.3 


+13:05:34 


2 


704 


27U 


28.54 


H[m0A[3H0S[323R[4i S EH 


NVSS J094925-200518 


09:49:25.8 


-20:05:18 


2 


706 


11 


27.25+ 




6C**1102+4329 


11:05:43.1 


+43:13:25 


2 


734 




28.17 




TXS 1545-234 


15:48:17.5 


-23:37:02 


2 


754 


27J1 


28.74 


OA. 3^3 OS[32l R[4l 


NVSS J152435-352623 


15:24:35.4 


-35:26:22 


2 


760 


Hp 


28.48 


TN J0920-0712 


09:20:22.3 


-07:12:18 


2 


762 


z!L 


28.48 




MG 1019+0534 


10:19:33.3 


+05:34:35 


2 


765 


298] 


28.57 


M| REH S[280l 


NVSS J111921-363139 


11:19:21.8 


-36:31:39 


2 


769 


11 


29.39 




TXS 0417-181 


04:19:43.6 


-18:01:56 


2 


770 


27JI 


28.64 


R[4lOS[320ll 


NVSS J213637-340318 


21:36:38.4 


-34:03:21 


2 


770 


38T 


28.39 




WNR J1338+3532 


13:38:15.1 


+35:32:04 


2 


772 


72] 


28.42 


R[24l 


7C 1807+6719 


18:07:13.4 


+67:19:42 


2 


78 


184] 


27.99 


4C +44.02 


00:36:53.5 


+44:43:21 


2 


790 


Bj 


28.98 




NVSS J140223-363539 


14:02:23.6 


-36:35:42 


2 


804 


up 


28.34 




WNH 1821+6117 


18:22:18.8 


+61:18:37 


2 


81 


2511 


27.94 




4C +03.16 


09:00:14.3 


+03:24:08 


2 


814 


2711 


28.48 




6C**0824+5344 


08:27:58.9 


+53:34:15 


2 


824 


66] 


28.27 




NVSS J094552-201441 


09:45:51.0 


-20:14:47 


2 


829 


n 


27. 12 1 ' 




NVSS J095357-203652 


09:53:57.4 


-20:36:51 


2 


829 


35] 


27.18+ 





Radio luminosity calculated assuming a spectral index a = -1.1. 



Table 3, cont 



Name 


RA(J2000) 


DEC(,T2000) 


z 1 


W.NhL 1717+6827 


17 


16 


56 


4 


+68 


24 


02 


2 


84 




01 


48 


28 


9 


+ 10 


28 


21 


2 


845 


T\/TT?f~' nn^9 9/11 


00 


54 


20 


8 


—23 


51 


30 


2 


86 


do i i go 1 07 
JjZ llOZ^O* 


11 


35 


05 





+37 


08 


41 


2 


88 


/IP 1 9/1 98 


13 


48 


14 


8 


+24 


15 


50 


2 


889 


TVQ ni^/1_L9^1 
_L _A_D UlO4i-|-Z01 


01 


37 


07 





+25 


21 


19 


2 


897 


IN V 00 JZloZoo-ooOolo 


21 


32 


39 





—33 


53 


19 


2 


900 


AC 98 ^8 


23 


51 


50 


1 


+29 


10 


29 


2 


905 


l\/TDp flQ/1 ^ 9/19 
IvirvO uy^o- Z4:Z 


09 


45 


32 


8 


— 24 


28 


50 


2 


922 




07 


47 


43 


7 


+46 


18 


58 


2 


926 


NVSS T1 51 090-352803 


1 ^ 
10 


1 n 

J_U 


911 

zu 





— 00 


9 S3 


Uo 


2 


yo / 


TYQ nQ^9 917 


09 


54 


20 


5 


—21 


56 


53 


2 


950 


"\^/"l\T"H" 17^8-1-^891 


17 


58 


48 





+58 


21 


41 


2 


96 


cp* nn9n_L/i/in 


00 


23 


08 


5 


+44 


48 


23 


2 


988 


WIN JU / 4 / -f-ODO-i 


07 


47 


20 


4 


+36 


54 


38 


2 


992 


DO Ui4!:0-r"0U±y 


07 


58 


06 





+50 


1 1 


03 


2 


996 


WIN JU/Ol-hoDUU 


02 


31 


11 


5 


+36 


00 


27 


3 


080 


"\\/"l\T T1 IT^^-L 1^/1 9/1 


10 


53 


36 


3 


+54 


24 


42 


3 


083 


do .izooUi-oyii 


23 


30 


24 





+39 


27 


1 1 


3 


086 


TNT T 1 1 1 9 9Q48 


1 1 


12 


2 





— 29 




06 


3 


09 


w Voo Juyo f ol-zloozl 


09 


/ 


51 


3 


— 21 


33 


2 1 


3 


126 


MRP 0^1 (S 9^7 


03 


18 


12 


1 


-25 


35 


10 


3 


142 


WN J0617+5012 


06 


17 


39 


4 


+50 


12 


55 


3 


153 


MRC 0251-273 


02 


53 


16 


7 


-27 


00 


10 


3 


16 


WN J1123+3141 


11 


23 


55 


9 


+31 


41 


26 


3 


216 


NVSS J230123-364656 


23 


01 


23 


5 


-36 


46 


56 


3 


220 


6CE 1232+3942 


12 


35 


04 


8 


+39 


25 


39 


3 


221 


WNH 1702+6042 


17 


03 


36 


2 


+60 


38 


52 


3 


223 


NVSS J232100-360223 


23 


21 


Of) 


6 


-36 


02 


25 


3 


320 


6C**0832+5443 


08 


36 


00 





+54 


33 


26 


3 


341 


B2 0902+34 


09 


05 


30 


1 


+34 


07 


57 


3 


395 


NVSS J094724-210505 


09 


47 


24 


5 


-21 


05 


06 


3 


377 


NVSS J095438-210425 


09 


54 


38 


4 


-21 


04 


25 


3 


431 



t 



Radio luminosity calculated assuming a spectral index a = -1.1. 
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Name 


RA(J2000) 


DEC(J2000) 


z 


NVSS ,7231402-372925 


23 


14 


02 


4 


-37 


2!) 


27 


3 


450 


TN J0205+2242 


02 


05 


10 


7 


+22 


42 


50 


3 


506 


TN J0121 + 1320 


01 


21 


42 




^1 3 


20 


58 


3 


517 


6C 1909+72 


19 


08 


23 




4-79 


20 


12 


3 


537 


4C 1243+036 


12 


45 


38 


4 


+03 


23 


20 


3 


560 


WN J1911+6342 


19 


I X 


49 


5 


+63 


42 


10 


3 


590 


MG 2141+192 


21 


44 


07 


5 


+ 19 


29 


15 


3 


592 


6C* 0032+412 


00 


34 


53 


1 


+41 


31 


32 


3 


658 


TN J1049-1258 


10 


49 


06 


2 


-12 


58 


19 


3 


697 


4C 60.07 


05 


12 


54 


8 


+60 


30 


51 


3 


791 


4C 41.17 


06 


50 


52 


2 


+41 


30 


31 


3 


792 


TN J2007-1316 


20 


07 


53 


2 


-13 


16 


44 


3 


837 


NVSS J231727-352606 


23 


17 


27 


4 


-35 


26 


07 


3 


874 


NVSS J021308-322338 


02 


13 


08 





-32 


23 


40 


3 


976 


TN J1338-1942 


13 


38 


26 


1 


-19 


42 


30 


4 


105 


TN J1123-2154 


11 


23 


10 


2 


-21 


54 


05 


4 


109 


8C 1435+63 


14 


36 


39 





+63 


19 


04 


4 


251 


6C 0140+326 


01 


43 


43 


8 


+32 


53 


49 


4 


413 


RC J0311+0507 


03 


11 


48 





+05 


08 


03 


4 


514 


TN J0924-2201 


09 


24 


19 


9 


-22 


01 


42 


5 


197 



Table 3, continued 



Log(L 50 o) 



28.83 
28.46 
28.49 
29.12 
29.23 

28.14 
29.08 
28.75 
28.94 
29.20 
29.18 

29.13 
28.71 
28.40 
28.70 

28.45 
29.40 
28.73 

29.49 
29.51 
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Radio luminosity calculated assuming a spectral index a = -1.1. 
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Table 4. HzRGs in order of right ascencion 



Name 






Name 






Name 






JVl Kj VJUUl-f-Uo4tO 


— 


~5kt\ — 




— — 


— — 




— — 


TTqq — 


T}~w<2 nmi niQ 

Jr ±V o U U ± l-Uzo 




nan 


1 _\. o U^i^ - 1 , ' . > 










„ 


MRC 0015-229 




010 


g^-i* * 0832 + 5443 




341 


4C -00 54 




360 


6C* 0020+440 




988 


60**0834+4129 




442 


PKS 1425-148 




or r 


NVSS J 002402-325253 




043 


LBDS 65W^172 




551 


8C 1435 + 63 




251 


6C* 0024 1 356 




161 


4C 12 32 






TXS 1436+157 




r 38 


4C 40 02 




428 


60**0854 + 3500 






NVSS J 144932-385657 




152 


MRC 0030-219 




168 


4C 40 22 






NVSS .1151020-352803 




937 


6C* 0032 + 412 




658 


4C +03 16 




814 


NVSS J 152435-352623 




760 


4C -4-44 02 




790 


6C 0901 + 54 






8C 1536 + 620 






WN J0040+3857 




606 


B2 0902 + 34 






TXS 1545-234 




7^4 


MRC 0052-241 




860 


TN J 0920-0712 






TXS 1 558-003 




r 20 


6C* 01064-397 




284 


TN J0924-2201 




197 


TXS 1647+100 






6C* 0112 + 372 


2 


535 


6C 0930 + 389 


2 


395 


WNH 1702 + 6042 




223 


6C* 0115 + 394 


2 


241 


MG 0936 + 0503 


2 


306 


TXS 1707+105 


2 


349 


6C* 0118 + 486 


2 


350 


MRC 0943-242 


^ 


922 


LBDS 53W002 


2 


390 


T N ,I0121-|-1320 




517 


NVSS J 0945 5 2- 2 01441 






W" ]-j 1717 + 6827 






6C* 0122 + 426 


" 




NVSS J 094724- 2 10505 






W NH 1736 + 6504 


~ 




BRL 0128-264 




348 


NVSS J094925-200518 


2 


706 


7C 1736 + 650 




400 


TXS 0134 + 251 


2 


897 


NVSS J094925-203724 


2 


427 


7C 1740 + 6640 


2 


100 


6C* 0135 + 313 


2 


199 


NVSS J 094949-2 13432 


2 


354 


MG 1747+1821 


2 


281 


MRC 0140-257 


2 


640 


NVSS J 09 5 22 6-2001 05 


2 


421 


MP J1755-6916 






6C 0140 + 326 


4 


413 


NVSS J 09535 7- 203652 


2 


829 


7C 1758 + 6719 


2 


700 


6C* 0142 + 427 


2 


225 


TXS 0952-217 


2 


950 


WNH 1758 + 5821 


2 


960 


MG 0148+1028 


2 


845 


NVSS J 09 543 8- 2 10425 




431 


MP J1758-6738 


2 


026 


6C* 0152 + 463 


2 


279 


NVSS J 09575 1-2 13321 




126 


7C 1802 + 6456 


2 


110 


MRC 0156-252 


2 


016 


NVSS J 101 008-383629 


2 


236 


7C 1807+6719 


2 


780 


TXS 0200 + 015 






MG 1019 + 0534 






4C 40 36 






TN J 0205 + 2242 






TN J1033-1339 




425 


W"NH 1821 + 6117 




810 


NVSS J 021308-322338 




976 


NVSS .J 103615-321659 




144 








TXS 0211-122 




340 


8C 1039 + 681 






^VN J1911 + 6342 






TXS 0214+183 




131 


gO* * 2045+4459 




571 


4C 48 48 






WN J0231 + 3600 




080 


TN J1049-1258 




697 


TN J 1941-1952 






MRC 0251-273 




160 


WN J1053 + 5424 




083 


TN J2007-1316 




837 


WN J0303+3733 


2 


505 


B2 1056+39 


2 


171 


MRC 2025-218 


2 


630 


RC J0311 + 0507 


4 


514 


TN J1102-1651 


2 


111 


TXS 2034+027 


2 


129 


MRC 0316-257 






6C**1102+4329 






MRC 2036-254 






MP ,J 0340- 650 7 


2 


289 


MRC 1106-258 




430 


MRC 2048-272 


2 


060 


MRC 0349-211 


2 


329 


6C 1106 + 380 


2 


290 


MRC 2052-253 


2 


600 


TXS 0355-037 


2 


153 


TN J1112-2948 




090 


MRC 2104-242 


2 


491 


MRC 0406-244 


2 


440 


WN J1115 + 5016 


2 


540 




2 


479 


TXS 0417-181 


2 


770 


TXS 1113-1 78 


2 


239 


NVSS J2 13238-3353 18 


2 


900 


TXS 0448 + 091 


2 


037 


4C 34.34 


2 


400 


NVSS .1213510-333703 


2 




TN J 0452- 1 737 


2 


260 


NVSS J 1 1 192 1-363 139 


2 


769 


NVSS ,1213637-340318 


2 


770 


4C 60.07 


3 


791 


3C 257 


2 


474 


MRC 2139-292 


2 


550 


PKS 0529-549 


2 


575 


TN J1123-2154 


4 


109 


MG 2141 + 192 


3 


592 


WN J0617+5012 


3 


153 


WN J1123+3141 


3 


216 


TXS 2202+128 


2 


704 


6C 0629+53 


2 


246 


B2 1132 + 37 


2 


880 


TNR 2254 + 1857 


2 


154 


4C 41.17 


3 


792 


6C 1134 + 369 


2 


125 


NVSS J230123-364656 


3 


220 


B3 0731 + 438 


2 


429 


MRC 1138-262 


2 


156 


4C +40.49 


2 


140 


PKS 0742 + 10 


2 


624 


B2 1159+39A 


2 


370 


MG 2308 + 0336 


2 


457 


WN J0747+3654 


2 


992 


B3 1204+401 


2 


066 


NVSS J231144-362215 


2 


531 


B3 0744+464 


2 


926 


4C 26.38 


2 


609 


NVSS J231402-372925 


3 


450 


6C**0746+5445 


2 


156 


6C 1232 + 3942 


3 


221 


NVSS J231727-352606 


3 


874 


TXS 0748 + 134 


2 


419 


WN J1242 + 3915 


2 


131 


NVSS J232100-360223 


3 


320 


6C**0743+5019 


2 


996 


4C 1243+036 


3 


560 


TXS 2319+223 


2 


554 


4C +39.21 


2 


119 


MG 1251 + 1104 


2 


322 


NVSS J232651-370909 


2 


357 


5C 7.10 


2 


185 


MRC 1324-262 


2 


280 


B3 J2330 + 3927 


3 


086 


5C 7.15 


2 


433 


WNR J1338+3532 


2 


772 


NVSS J233034-330009 


2 


675 


B3 0811 + 391 


2 


621 


TN J1338-1942 


4 


105 


TXS 2332+154 


2 


481 


5C 7.223 


2 


092 


4C 24.28 


2 


889 


MG 2348+0507 


2 


019 


6C**0824+5344 


2 


824 


4C +43.31 


2 


149 


4C 28.58 


2 


905 


5C 7.269 


2 


218 


PKS 1357+007 


2 


671 


TXS 2353-003 


2 


592 
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